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^•^ Abstract 

We use chiral perturbation theory (ChPT) to calculate the ji^d scattering length with an accu- 
racy of a few percent, including isospin-violating corrections in both the two- and three-body 
sectors. In particular, we provide the technical details of a recent letter [111, where we used 
data on pionic deuterium and pionic hydrogen atoms to extract the isoscalar and isovector pion- 
^ ■ nucleon scattering lengths a^ and a . We study isospin-breaking contributions to the three-body 

P5 \ part of a,j-d due to mass differences, isospin violation in the jiN scattering lengths, and virtual 

photons. This last class of effects is ostensibly infrared enhanced due to the smallness of the 
»^ \ deuteron binding energy. However, we show that the leading virtual-photon effects that might 

^ ■ undergo such enhancement cancel, and hence the standard ChPT counting provides a reliable 

0^ \ estimate of isospin violation in ajj-d due to virtual photons. Finally, we discuss the validity of 

^^ ■ the Goldberger-Miyazawa-Oehme sum rule in the presence of isospin violation, and use it to 

determine the charged-pion-nucleon coupling constant. 
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s_j ■ 1. Introduction 



Hadron-hadron scattering lengths are fundamental quantities characterizingthe strong inter- 
action, and are slowly becoming accessible to ab initio calculations in QCD ||2l|3|]. Among them, 
of particular interest are pion-hadron scattering lengths: the chiral symmetry of QCD and the 
Goldstone-boson nature of the pions dictate that they are small I4l|, while their non-vanishing 
size is linked to fundamental quantities like the light quark masses and condensates. For exam- 
ple, the combination of two-loop chiral perturbation theory (ChPT) and Roy equations resulted 
in very precise predictions for the pion-pion scattering lengths ||5|] 

flo = (0.220 ±0.005)M^', «„ = (-0.0444 + 0.0010)M^\ (1.1) 
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which were essential to confirm the role of the quark condensate as the leading order parameter 
in the spontaneous breaking of chiral symmetry j^ . 

In the case of pion-nucleon scattering, chiral symmetry predicts that the isoscalar scattering 
length fl^ is suppressed compared to its isovector counterpart a . In particular, the low-energy 
theorem for a ||4|,|7|] 

+ O(mI)^SQ-10-^M„^ (1.2) 



87r(l + MJm^)Fl 

receives corrections only at third order in the pion mass and its prediction is numerically very 
close to the full result. Meanwhile, the expansion of the isoscalar scattering length 171] 

with the pion decay constant F-„, the axial charge of the nucleon ^a, and low-energy constants 
(LECs) c,, stands in marked contrast: the leading order vanishes — leaving a^ as a measure of 
the explicit breaking of chiral symmetry — and at sub-leading orders poorly determined LECs 
and huge cancellations between individual terms limit the predictive power of the expansion. 
Experimentally, lack of :7r" beams and neutron targets makes direct pion-nucleon scattering ex- 
periments impossible in some charge channels, complicating a measurement of a^. In the isospin 
limit the tt"/? scattering length is purely isoscalar, and corrections to the isospin limit are well- 
controlled for this quantity. The best hope for access to a+ in the jiN sector therefore lies in 
precision measurements of threshold neutral-pion photoproduction ISllst]. But, until the advent 
of such measurements, extractions of a^ from jiN scattering data suffer from large uncertain- 
ties. Different phase-shift analyses yield values covering a wide range from -10 ■ lO^^M^' to 
+5 ■ lO^^M^' llOfl. Indeed, the combination of data and theory has, until now, lacked sufficient 
accuracy to even establish definitively that a^ + 0. 

A precise determination of jiN scattering lengths improves our knowledge in many areas; 
two particularly important examples of this are the following. First, a^ is one of several inputs 
to dispersive analyses of the pion-nucleon cr-term lllil . which measures the explicit chiral sym- 
metry breaking in the nucleon mass due to up and down quark masses, and is, in turn, connected 
to the strangeness content of the nucleon. Second, a~ serves as a vital input to a determina- 
tion of the pion-nucleon coupling constant via the Goldberger-Miyazawa-Oehme (GMO) sum 



rule Ill2ll . While the uncertainty in a' is much smaller than that in a^, it still contributes signif- 
icantly to the overall error bar on the sum-rule evaluation 1 13111411 . This latter example is thus 
one of several where data on pion-nucleon scattering affects more complicated systems like the 
nucleon-nucleon (A^A^) interaction, and hence has an impact on nuclear physics. 

In view of the difficulties concerning both direct experimental access and the convergence of 
its chiral expansion ( 11.31 ). data on hadronic atoms have become the primary source of information 
on a^ fll5n . In these systems, the strong interaction modifies the spectrum compared to pure QED 
by shifting the energy levels and introducing a finite width to the states. Both effects are sensitive 
to threshold pion-nucleon scattering. In this way, new information on pion-nucleon scattering 
lengths has become available due to recent high-accuracy measurements of pionic hydrogen {nH) 
and pionic deuterium {nD). In the case ofnH, the latest experimental results 11611 are 

ei, = (-7.120 ±0.012) eV, r,, = (0.823 + 0.019) eV, (1.4) 

for the (attractive) shift of the 1 s level of ttH due to strong interactions and its width. The shift 
of the ground state is related to the n^p scattering length a„-p, while the width gives access to 



the charge-exchange scattering lengtha^™ = a^-p^^^o,, 11511 . More precisely, £is is related to a„-p 



through an improved Deser formula 11711 



ei, = -2aX«--p(l + K. + S^), (1-5) 

where a - e^jAji, ^h is the reduced mass of jiH, K^ = 2a(l - \oga)^Ha7rp, and iS™*^ - 
26^h(0)/^h(0) = 0.48% is the effect of vacuurn polarization on the wave function at the ori- 
gin ilST . The width, in turn, determines a'^^'^ via 11911 



Tu = 4a'f/„pi{l + ;^)«-P'(l +Kr + 5^1 (1.6) 

with/fr = 4o!(^ -iog a)fiHajT- p+2fiH{inp+ M„ -m„- M„o)(a„o„)^ , nip, m„, M„, andM;^o the masses 
of proton, neutron, charged and neutral pions, respectively, p\ the momentum of the outgoing 
nn^ pair, and the Panofsky ratio 120(1 

P=— -^ ^ = 1.546 + 0.009. (1.7) 

cr(;r p — * ny) 

Similarly, the (repulsive) strong shift e'^^ of the 1 s level of nD yields the real part of the tt li 



scattering length Re a„-d via 112 111 



ef, = -2aVoRe a„-d(l + K,o + 6J), (1.8) 

where fip is the reduced mass of nD, K^d = 2a{ 1 -log Q')yU/)Re a„-d, and 5™*^ = 25'P£)(0)/*I'z)(0) = 
0.51% yj]. 

In the isospin limit, the level shift of ttH is sensitive to a^ + a , whereas the width is solely 
determined by a . In this way, data from ttH alone permit, in principle, an extraction of the 
ttN scattering lengths. However, the chiral suppression of a^ makes it very sensitive to isospin- 
violating corrections (see Sect. |2|, such that additional experimental information — in particular 
from isoscalar nuclei as they provide better access to a^ — are essential to check the systematics 
and potentially improve the accuracy of the scattering-length determination. To this end, we split 
the jT^d scattering length into its two-(7rA^) and three-(7rA^A^) body contributions 



Rea„-d = a%+a% (1.9) 



where the former is related to a^ via 



2^P.~. 



«f-'/=T^(«'+A«')- (1-10) 

Here the difference between a^ and a^ as well as Afi^ are determined by isospin-violating cor- 
rections (Sect. 121) and 

nip niA 

with the deuteron mass m^. Once isospin breaking in the two-body sector is under control, 
we therefore have to develop a theoretical description of ajj that finally allows one to exploit 
information on jiD at the same level of accuracy as in nH, which requires that we can compute 
'^TTd '■° ^^ accuracy of better than 10%. As we shall discuss in this work, this proves to be 



possible, and a combined analysis of the data ( 11.41 ) on nH and the recently remeasured level shift 



irniDUM 



efj = (2.356 + 0.03 l)eV (1.12) 

then yields the determination of a* and a' of unprecedented accuracy in yj]. Qlie width of nD 
is governed by n^d — » nn (BR - 73.9 %) and n^d — > nny (BR = 26.1 %) 12311 . such that no 
additional information on threshold nN physics is provided.) The main purpose of this paper is 
to provide the details of the calculation of the three-body part of ajj-j, which we decompose as 

where a'^'^p^'^ involves two-nucleon or A-isobar intermediate states, a^^ represents virtual-photon 
corrections, and cf^^ denotes "strong" diagrams, i.e. essentially all other contributions in the chiral 
expansion (the definition of each class of diagrams can be found in Sects. |4]-[6|. 

The paper is organized as follows: we first briefly review isospin-violating corrections to the 
TiN scattering lengths in Sect. |2l Then, we summarize the hierarchy of diagrams contributing to 
a^Ij in both the isospin-conserving and the isospin-violating sector in Sect. [3j before discussing 
strong, virtual-photon, and dispersive + A contributions in detail in Sects. |4]|5] and|6l A reader 
not interested in the details of the calculation may skip Sects. |4}{6]and proceed to Sect.|7] where 
we summarize our main conclusions concerning three-body contributions to the ji^d scattering 
length. The consequences for the nN scattering lengths and the jiNN coupling constant are 
presented in Sects. [8]and|9] We conclude in Sect. [TO] Various details of the calculation are 
provided in the appendices. 

2. Isospin violation in the nN scattering lengths 

Before turning to the calculation of a^_^j, we review isospin-violating corrections to the nN 
scattering lengths, which provide an essential input to the present analysis. The scattering lengths 
in the isospin limit for all eight channels can be written in terms of a^ and a~ as 

^n~p = ^7:~ p—*n~ p ~ ^n'^n = ^n^n—^n'^n —CI -r Q , 



ccx ccx 



■^TT^p—^n^'p ^n'n — ^n^'ii—^n^'ii 



■yl2c 



(2.1) 



To extract a^ and a from hadronic-atom data, we need to relate the scattering lengths in partic- 
ular charge channels to those in the isospin limit, i.e. we need the corrections 

^a„-p = a„-p - (fl+ + fl"), ^a„-„ = a„-„ - (fl+ - a"), Aa^"^ = 0"^:^^ + yfla-. (2.2) 

These corrections are generated by the quark mass difference m^ - m^ and electromagnetic in- 
teractions. They can be calculated systematically in ChPT, and have been worked out at next-to- 



leading order (NLO) in the chiral expansion in Il24l - j26ll 



In those works, and throughout this study, the counting m^ - m^ ~ e^ is used, i.e. electro- 
magnetic and quark-mass effects are assumed to contribute at the same order. This counting is 
phenomenologically rather successful. The prime example is the nucleon mass difference, to 



which — according to the evaluation of the Cottingham sum rule in |27!l — the quark mass dif- 
ference and electromagnetic interactions contribute (2.1 + 0.3) MeV and (-0.8 + 0.3) Me V, re- 



spectively. (This result is consistent with recent determinations from the lattice Il28ll29ll and from 
charge symmetry breaking in pn -^ dn^ 13011 .) A similar picture emerges from the kaon mass dif- 
ference, where — depending on the assumptions about violation of Dashen's theorem i3ll433n — 



quark-mass effects are a factor 2-3 larger than electromagnetic ones. It is also instructive to look 



at tree-level contributions to isospin violation in nN scattering 112411 : a^j-p - aj:+„ and aj^+p - a„-„ 
are purely electromagnetic, a,fip — a„o„ is solely due to nid - m^, while both effects are of the same 
.. _ g^™. Similar conclusions can be drawn from tree-level isospin breaking in the nK 



scattering lengths 113411 . where the corrections for some channels are purely electromagnetic, for 
some purely quark-mass induced, and for some due to both effects, sometimes the former being 
a factor of 2 larger, sometimes the latter. 



First of all, the major consequence of the leading-order (LO) isospin breaking in ChPT 1I35I1 

is that it is impossible to directly extract a^ from hadronic atoms. Only the combination 

is accessible, and a^ itself cannot be obtained absent input on the LECs ci and /i from other 



sources (the full list of LECs relevant for the present work is given in Appendix A i. If the stan- 
dard single-nucleon-sector counting e ~ p is employed, then these isospin-violating effects are 
actually of the same size as the piece of 5^ that would be present in the isospin limit, ci en- 
ters these effects because its contribution to a^ is proportional to M\, and /i features in the 
electromagnetic contributions to nij, and m„. Estimates of these constants will be discussed in 
Sect. [81 

Since only a^ can directly be extracted, it is convenient to work with 

Aa^p = a^p - (5+ + a"), Aa,r„ = ajm - (5^ - a') (2.5) 

instead of Aa„-p and Aa„-„. The results relevant in the present context may then be written as 

Aa„-p — Aa^ + Aa^, Aa„-„ — Aa^ - Aa", 



^« = —TT- - T-r TTTT^ 3 + log — f -H — 5-c/^ + -—-^ 1 -H 41og2 + 3 log -f 






M„h„ 



j(-'«i^<-'«i)-f)-3^(— f 



1927:2^4 



Fl 



4 + ^ [pls'i - 2^3 + n + f (M + 4)) \ (2.6) 



where 



An = mn - mp. (2.7) 

The apparent dependence on the renormahzation scale yu is canceled by the scale dependence of 



the LECs, whose definition is briefly reviewed in Appendix A (for more details we refer to 112611 V 



Estimating the LECs as in 02411 yields 



Aa+ = (-3.3 + 0.3)- \^-'^M-\ Aa- = (1.4 + 1.3) ■ lO^^M^', 
Afl;,-;, = (-2.0+1.3)- IO^^^M;', Aa^^-'), = (0.4 + 0.9)-10"^^M;'. (2.8) 

In principle, one could also define an 5" in which some of the LECs appearing in Aa" and 
AaJ^':" in the same way can be absorbed (namely /2, (ig, A:p and k^), similarly to the definition 
of a with respect to ci and /i . In this way, the constraints of kH and nD on 5^ and 5" would 
be considered and only in the end the estimates for the LECs inserted. The advantage of this 
alternative procedure is that the dependence on the LECs is more transparent and correlations 
between the three constraints better under control. However, we have checked that the results 
obtained in such an approach differ only marginally from what we present here. In order to keep 
the discussion as simple as possible we work in terms of 5^, A5^, a , and Aa . 

3. Hierarchy of three-body operators and Weinberg power counting 

3.1. Isospin-conserving operators 

So far no counting scheme is known that permits consistent, realistic, and simultaneous con- 
sideration of the two- and three-body operators that contribute to n^d scattering. For example. 



in the original power counting by Weinberg 0361 13711 the leading two-body operator (propor- 



tional to a^) appears formally at one order lower than the leading three-body terms shown in 



Table [T] However, it has been known for years (see e.g. 0131 14411 and references therein) that 
the double-scattering diagram (the first diagram in the first row of Table [T]i alone is close to the 
experimental value of the scattering length, whereas the term proportional to a^ is significantly 
smaller This drawback of the counting scheme may bring into question the theoretical uncer- 
tainty estimate obtained within ChPT and thus the rehability of a ChPT extraction of the nN 
scattering lengths from pionic deuterium. On the other hand, practical calculations demonstrate 
that Weinberg's power counting still works quite well once it is applied to two- and three-body 
operators independently — in spite of the difficulties that this power counting has in accounting 
for the relative size of these two classes of contribution. In particular, it was shown in fi% that an 
application of the Weinberg scheme allows one to systematically account for strong three-body 
contributions to aj^d to very high accuracy. In this work we demonstrate (see Sect. |5]l that Wein- 
berg power counting is also fully in line with the actual size of the isospin-violating three-body 
contributions. Since isospin breaking in the two-body sector is also well under control 02411 . this 
permits a precise extraction of 5^. 

Therefore, in what follows we consider a power counting within the class of three-body con- 
tributions. All such diagrams are ordered using the Weinberg scheme, with their order quoted as 



Chiral order 



Three-body operator 



Reference 



LO = 0(1) 



> 



[36,37] 



NLO = 0(p) 



mi 



Effect of nucleon recoil in LO diagrams 



N^/^LO = 0{p'/^) 



t 



[42] 



Effect of nucleon recoil in LO diagrams 



[43J 
[40,41J 



N-LO = Oip^) 




+ 



Table 1 : Hierarchy of isospin-conserving three-body operators within Weinberg power counting. Solid (open) circles 
correspond to leading (sub-leading) vertices, grey blobs indicate the deuteron wave functions, and the black ellipse 
coiTesponds to NN interactions in the intermediate state. Solid single, solid double, and dashed lines correspond to 
nucleons, A(1232)-isobars, andpions, respectively. 



the predicted size in that counting relative to the leading, double-scattering term. The isospin- 
conserving three-body diagrams are illustrated in Table[T] The table also quantifies their relative 
contribution according to the small expansion parameter p. We note that the ordering in Table [1] 
goes beyond naive dimensional analysis, since known enhancements and contributions at frac- 
tional orders in the expansion parameter are already taken into account (see below). The goal 
is to include all three-body operators up to one order lower than the contribution of the lead- 
ing unknown {N'^N)-nn contact term, which appears at next-to-next-to-leading order (N-LO). 
Its contribution cannot easily be determined from data, and is a key source of uncertainty in our 
result. Given that 0{p) ~ x - M„/mp, we anticipate an accuracy of a few percent for threshold 
TT^d scattering. 

This expectation is substantiated by the sensitivity of our integrals to the choice of the 
deuteron wave function (see Sect. 14.4b . Convolving the operators of Table [T] with different wave 
functions derived from chiral and phenomenological A^A^ potentials we find a variation in the 
results of about 5 %: an independent estimate of the contact term's effect. The explicit cutoff 
dependence of the individual diagrams at LO and NLO has been extensively discussed in the 
literature 1391 I45l - l47ll and it is now well established that, for deuteron wave functions based on 



the one-pion-exchange interaction, the results become cutoff-independent in the limit of a large 
cutoff. 

The hierarchy of three-body operators contributing up to N'LO is shown in Table [1] In 
addition to the double-scattering diagram there are also LO diagrams involving 3nN^N and 4n 
vertices, which individually depend on the parametrization of the pion fields, while only the sum 



is parametrization independent. The effect of these diagrams is numerically irrelevant Il36ll38ll . 
The reason for that was understood in ll39ll to be due to an accidental suppression of the spin- 
isospin matrix element, which appears to be more than one order of magnitude smaller in these 
diagrams than that in double scattering, although the momenta in the diagrams are in line with 
Weinberg power counting. 

The o pera tors at NLO involve sub-leading vertices and were shown to cancel amongst them- 
selves in 0811 . In addition, at NLO there is a triple-scattering term that requires some care. The 
problem is that the actual size of this diagram is enhanced as compared to the estimate based on 
dimensional analysis, which predicts that this diagram contributes only at N^LO. Specifically, 
it was shown in |i390 that the long-ranged (infrared) part of this diagram is enhanced numeri- 
cally by a factor of 7t~ whereas the rest behaves in accord with Weinberg counting. The origin 
of this enhancement was associated in ll39ll with the special topology of the diagram consisting 
of two consecutive pion exchanges with Coulombic-type pion propagators. It is interesting to 
note that enhancements by factors of n were already observed to emerge also in pion-loop con- 
tributions to the A^A^ potential ll48ll . the scalar nucleon form factor |49ll. the jfi photoproduction 



amplitude 115011 . and even isospin violation in ttN scattering 1124 12611 itself, cf. (I2.6l l, from sim- 
ilar topologies as those discussed here. A deeper understanding of when these dimensionless 
factors appear would be very desirable. For tt c/ scattering this sort of numerical (not para- 
metric) enhancement may also appear in topologically analogous diagrams that belong to the 
so-called multiple-scattering series. This leads to concerns regarding the quadruple-scattering 
term: although it formally appears only at N'^LO, which is far beyond the edge of the theoretical 
accuracy, a potential numerical enhancement needs to be carefully studied since it may affect the 
uncertainty estimate. In fact, the whole class of multiple-scattering diagrams can be summed up 
to all orders, and we find the effect from quadruple-scattering and higher diagrams to be negli- 
gible, see Sect. l4.3l for more details. Thus, this class of potentially dangerous diagrams does not 
affect our uncertainty estimate. 



In addition, starting from NLO, nucleon-recoil effects to the leading double-scattering oper- 
ator have to be taken into account. The nucleon recoil has been widely studied in the literature 
both phenomenologically, see e.g. i5lll52ll . and within effective field theory (EFT) 140114 Up 311 . 



In 114 111 it was demonstrated how one can calculate the recoil effect to all orders in a systematic 
expansion within EFT. At leading order nucleons are considered as infinitely heavy, such that 
the pion is scattered off two fixed centers. At NLO the nucleon kinetic energies enter and the 
static pion propagator needs to be replaced by the full propagator corresponding to the three- 
body ;rA^A^ intermediate state. In the regime where all momenta in the diagram are of order of 
Mj^ the nucleon-recoil effect is purely perturbative, and thus it can be calculated by expanding 
the nucleon kinetic energies using standard heavy-baryon techniques. In this regime the recoil 
effect contributes at integer powers in the expansion, i.e. at 0{p), O(p^), etc., relative to the lead- 
ing, double-scattering effect. In addition, there is also a non-perturbative regime in which the 
three-body propagator goes to zero — the regime of the three-body singularity. In this regime the 
pion kinetic energy is of the order of the nucleon recoil, such that pion momenta appear to be 
suppressed by ^M„lm^ [^ compared to the characteristic momenta in the deuteron. Thus, the 
expansion of the double-scattering diagram contains half -integer powers of M^/nip due to the 
recoil effect. Note that the potentially largest isovector recoil correction at order 0{p^^^) fully 
determined by the small scales vanishes exactly as a consequence of the Pauli principle, which 
prohibits the A^A^ intermediate state to be in an S -wave in this case i40ll4I[l . Furthermore, it was 



demonstrated in 1(41(1 that the recoil effect for n'd scattering is relevant only at orders 0(p) and 
0{p^^^), whereas the contributions at higher orders are already negligible. This can be considered 
further support for the Weinberg counting scheme and our resulting uncertainty estimate. 

In order to achieve the desired accuracy we must, in addition, compute the dispersive cor- 
rections as well as the contributions that emerge due to the expUcit treatment of the A(I232) 



resonance ll42ll43ll . which enter at order 0{p^^^), see Sect. |6]for more details. Both effects in- 
volve new scales. The dispersive corrections due to the process nd — > A^A^ — > nd are linked 
closely to pion production in A^A^ collisions, where, as a consequence of the large momentum be- 
tween the A^A^ pair, the expansion parameter is ^^JM^Jm^. The A-resonance in the A-less theory is 
hidden in the low-energy constants c, of nN scattering and it contributes to tt li scattering through 
the so-called boost (Fermi motion) correction 1381 39ll rl This two-body correction was shown 



to be significant but strongly model-dependent 113811 . although in a more refined treatment 1I39I1 
the model dependence was shown to be smaller In any case, the A-nucleon mass difference is 
just about twice as large as the pion mass and it is profitable to include the A(1232) dynamically 
and in this way increase the breakdown scale of the theory. When this is done the value of the 
relevant LEC C2 is reduced by almost an order of magnitude. Therefore the inclusion of the A as 
an expUcit degree of freedom allows one to reduce the model dependence and to achieve a faster 
converging series through the explicit calculation of a certain class of diagrams 11391 14311 . The 
residual boost correction is negligible in the A-full theory y. 



' Here and below we identify the nucleon mass with the mass of the proton, apart from instances where we explicitly 
discuss the impact of the nucleon mass difference An on our results. 

-Note that the LEC ci and the Hnear combination of LECSC2+C3 contribute to the ;rA' scattering length a* and through 
that also to n^d scattering. However, neither C[ nor C2 + c's are affected by the A-isobar up to order O(p^) 1541 , although 
the values of ci and C3 individually are strongly saturated by the A and thus change significantly when considering the A 
as an explicit degree of freedom I54ll5al . 



3.2. Isospin-violating operators 

Two-body isospin-violating corrections were already discussed in Sect. |2l Here we will 
follow the logic of the previous subsection and apply the Weinberg power counting to the 
isospin-violating three-body corrections. Thus, we are going to discuss the hierarchy of isospin- 
violating three-body operators relative to each other, and their relative suppression compared 
to the leading-order isospin-conserving operators. This is sufficient to perform a high-accuracy 
calculation of the n^d scattering length. 

There are three classes of isospin-violating three-body contributions: first, isospin-breaking 
corrections appear in nNN propagators due to pion and nucleon mass differences; second, 
isospin-violating jiN interactions can occur in the diagrams introduced in Sect. 13.11 The opera- 
tors corresponding to these two classes of isospin-violating mechanisms are marked as crossed 
circles in Tabled Since isospin violation in hadron masses and nN interactions can occur due to 
both electromagnetic and quark-mass effects both these classes could be either ~ nii — m^ or ~ e^. 
The third class is a purely electromagnetic effect: a new set of diagrams involving (low-energy) 
virtual photons (see Tableland Sect.|5]for more details). 

At leading order in isospin violation diagrams of this third class that involve a virtual photon 
and one insertion of the isospin-conserving jiN vertex occur. These are represented by the first 
row of diagrams in Table |2] At the same order effects due to the pion mass difference in the 
leading-order isospin-conserving diagrams enter (the second row of diagrams in Table|2]i. These 
effects can be computed by working on the particle basis for the pion intermediate states in 
the leading-order three-body diagrams (see first row of Table [1]), and explicitly keeping track 
of charged and neutral pion masses there. However, when this is done the double-scattering 
diagram must be treated in a special way, because its jiNN intermediate state, and associated 
three-body cut, means that the pion-mass-difference effect in this graph generates contributions 
not just at O(e^), but also at order (9(e^/7'^^), O(e^p), etc. The double-scattering graph with one 
insertion of the pion mass difference is therefore shown in Table |2] at LO, N'^^LO, and NLO. 
The piece of this graph which is LO in isospin violation is suppressed by e^F'^jM^ compared to 
the corresponding isospin-conserving diagrams at LO, as are the other diagrams listed in the first 
two rows of Table|2] Diagrams involving NNy intermediate states are also of this size (see third 



line of Table|2]i, but these are included in the calculation of the dispersive corrections in 14211 . and 
so will be accounted for in Sect. |6] 

Next in importance are the non-analytic terms which result from the inclusion of the pion 
mass difference in the nNN propagator of the double-scattering diagram. These yield the N'''^LO 
contribution of Table|2] in full analogy with the effect of nucleon recoil in the isospin-conserving 
case. 

The operators at NLO are suppressed by 0{e^p) compared to the three-body isospin- 
conserving operators at LO, and, given the smallness of the expansion parameter, are irrelevant 
for our present purposes. (Furthermore, for full consistency with the power counting, the inclu- 
sion of NLO corrections to the three-body isospin-violating operators would require the calcula- 
tion of N^LO isospin-violating two-body corrections, something that has not yet been achieved.) 
Therefore, for this study, it is necessary only to calculate all isospin-violating corrections to the 
n^d scattering length up to N'''^LO. 

However, due to the appearance of new scales in the three-body problem, there might be 
some higher-order operators that are enhanced, even though formally they only appear beyond 
0{e^p^^^). We identify and investigate these contributions explicitly: 

a) To account for all effects related to the three-body cut in the double-scattering diagram we 
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Chiral order 



Three-body operator 



Reference 



LO = 0(e^) 



'u- 



-V 



t 




[42J 



N'/^LO = 0(e'-p'/^) 



NLO = 0(e2p) 



+ 



N^LO = 0(e^p^) 



t 



+ 



Table 2: Hierarchy of isospin- violating three-body operators within Weinberg power counting. The suppression of these 
operators is given with respect to the isospin-conserving diagrams at LO (cf. Table0. Isospin violation appears either 
due to the inclusion of virtual photons or due to mass differences and electromagnetic effects marked by crossed circles. 
Note that diagrams with intermediate NNy states were already considered in (4^. The contribution of these diagrams is 
included in the result for the dispersive coiTections which we will adopt from 1431 . 
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contributions to n^d scattering. 
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keep all terms proportional to the pion and nucleon mass differences in the jiNN propagator 
unexpanded (see Sect. |4]i. In particular, we include the nucleon mass difference in the 
propagators to have the jfinn and n'pn thresholds at the proper places, although this is 
formally an NLO effect^ 

b) Due to the large size of the double-scattering diagram in the isospin-conserving case we 
include isospin violation in the nN vertices in this diagram. This effect also starts at NLO. 

c) We study certain virtual-photon corrections to the double-scattering process (formally ap- 
pearing at N^LO). The presence of virtual photons enhances the region of small momenta 
in these diagrams, such that the integrals become infrared divergent in the limit of van- 
ishing deuteron binding energy. The finite binding energy of the deuteron renders these 
diagrams finite, but the resulting contribution is potentially enhanced. In view of the fact 
that double scattering is numerically by far the dominant contribution to n^d scattering, 
these virtual-photon corrections could become relevant for our study. This effect is dis- 
cussed in detail in Sect.|5] 

Note that the isospin-violating three-body mechanisms up to N'^^LO are of purely electro- 
magnetic origin, while isospin violation due to the quark mass difference appears only in higher- 
order corrections, e.g. a) and b). Recall that isospin violation in the pion mass difference is 
predominantly an electromagnetic, not a quark-mass effect. There is no piece of this quantity 
~ BJd - nia at LO in isospin violation. Nucleon-mass-difference effects enter contribution a) only 
at O(e^p). In the case of b), the isospin-violating nN interactions which appear there include 
terms proportional to the quark mass difference — as reviewed in Sect. |2] But this whole class of 
diagrams involving isospin-breaking pion-nucleon vertices does not start until NLO. 

We demonstrate in the subsequent sections that it is indeed the case that the additional correc- 
tions a)-c) beyond N'^^LO are significantly smaller than the estimated theoretical uncertainty of 
the full analysis. The explicit calculation of these higher-order corrections provides an additional 
test of our counting scheme and uncertainty estimate. 

4. Strong contributions to n'd scattering 

4.1. Double scattering 

We start our discussion with the double-scattering diagrams {d\) and (A) (cf. Fig.[T]i. The 
diagram {d\) appears already at LO and contributes to all higher orders due to the effect of 



^At first order in isospin breaking we have An = -4Bc5(m(j-m„)+/2e^F^ and A;t = 2Ze^F^,withc5 ~ fi ~ l/wp and 
Z = 0(1). These quantities then enter the nNN propagator in the combination p = 2M;tAn - A;t, cf. 14. U . Assuming that 
the electromagnetic and quark-mass-induced contributions to An are of the same size, one finds M„A!^/A,r ~ ^^rr/'^p ~ 
0(p). Therefore the nucleon-mass-difference contribution top is suppressed by one chiral order compared to that coming 
from A;t- After accounting for the modification of the chiral counting due to the presence of the jtNN cut we find that A^ 
contiibutes to a"^"' at 0(e^p), whereas A;; affects the scattering length already at 0(e^). 
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nucleon recoil, as was discussed in the previous section. The diagram (d2), however, gives rise 
to a three-body contribution only if nucleon recoil is included. In the limit of static (infinitely 
heavy) nucleons the pion loops in diagram (dz) are already subsumed in the nN scattering lengths, 
since in this limit (dz) is nothing but an ordinary loop correction in the chiral expansion of 
the ttN scattering lengths. In this way, the contribution to n^d scattering from the part of (c/a) 
corresponding to static nucleons is always included in the two-body term a|^_^ proportional to a^, 
see (II.9I 1 and (II. 10b . Thus, to obtain an additional three-body correction we need to investigate 
the effect of embedding the nN amplitude into the nNN system. In this way we explicitly see 
one of the recoil effects; in order to treat the three-body dynamics properly we must replace 
the contribution of the two-body (nN) cut by that of the three-body (nNN) cut. For the isospin- 



conserving case the procedure for this was established in 1140114 ID . The goal of this section is to 
extend it to the isospin-violating case. 

The isospin-violating corrections occur in the diagrams of Fig.[T]through the different masses 
of particles in the intermediate states and through isospin-violating corrections to the nN scatter- 
ing lengths. Note that in the calculation of the double-scattering diagrams one can safely omit 
all isoscalar terms, since the term proportional to (5^)^ is tiny compared to (a^)^, while the term 
proportional to the combination a^a^ cancels. Therefore we calculate the diagrams of Fig. [T] 
keeping only the isovector nN scattering amplitude, retaining the isospin-violating correction 
Aa-. 

In this way, we can use the form of a'^"' that is correct in the presence of a pure isovector nN 
interaction. (The inclusion of recoil effects in the isoscalar case is discussed in [41].) Following 
the procedure described in |J40|], we obtain: 



aidOHd.J ^ ^static ^ ^static ^ ^cut ^ ^^(2)^ 



2/ 1 \ ^.static ^2 1 1 / ^' 



static _ 7;2/ \ ^static 



a = -a ( ^r ), a'^un = «^(— ( 



q2/' "^'LO \qzv^^+„,p 



A«^^) = «L [d\ (^ - -^^l (4.1) 

J Vq- + 5 q^+6 + p' 



where 



<5, 



Pl+P2\ ,, 3, ^ '^qq' 



/ PI + P2 \ ~ '^qQ" 
2cOp^-J,A£ + ^- , 5 = 5pp_q, 6^ , p = 2M;tAn - A;r, 

\ znip I nip 

^Ml+q\ </(q)) = Jd^Vd'q^'Xp - q)/(q)*I'(p), (4.2) 



P1.P2 



and e is the deuteron binding energy. The nN scattering lengths in ( 14.1b are defined as 

i:2 , -, . £2 



_L_((«- + A«-)2 + \ («- f ) = J^^{2{a-f + 2a- Aa- - V2a-Aa- +■■■), 



-a" 



^^ / x2 ^^ 



cex r, 2 



(<- ) = -2f((«"^' - ^«"^«r. + ■■■)' (4.3) 



In^^i^^P' n^U 



13 



where the elhpses contain higher orders in isospin violation. We use a normaUzation of the 
deuteron wave functions ^^(p) where 



/ 



d^p^P^Cpmp) = 1. (4.4) 



The individual terms in ( 14.1b can be interpreted as follows: a'^'^"'^ corresponds to (di) evaluated 
with a static pion propagator, and is numerically by far the dominant contribution. Recoil cor- 
rections to the static pion propagator are incorporated in fl^LO' while a'^"' comprises effects due 
to the three-body ;7r"«« and n^pn cuts. Aa*^' emerges as an isospin-violating correction in this 
rearrangement, which in the end does not constitute a true three-body effect — as indicated by the 
absence of the deuteron wave function. In the isospin limit p = A„ = A^ = Aa' - Aa'^^J^ - 



(14.11) reduces to the result derived in 1401 . 

Our power counting is based on dimensional analysis assuming all integrals scale only 
with M„. In fact, the integrals in ( 14.11 ) involve other scales too: '\/M„e, ^M„lm^ M„ — due 
to the three-body cut — , and ^m^e, thanks to the deuterium wave functions. The appearance of 
the first two scales becomes apparent if one realizes that in the regime of the cut the pion kinetic 
energy is of the order of the nucleon recoil so that pion momenta are suppressed in this regime. 

At first glance, the presence of a three-body cut in the integral for a'^"' makes it appear to be 



enhanced over its naive ChPT order by ^m^jM,, 15311 . Indeed, it was shown in 114 ill that the full 



result for the double-scattering diagrams can be expanded in half-integer powers of;^' - M„lm.p 
^WiH(d2) ^ ^static ^^1/2^^ ^^^^ +x"-a, + ■ ■ ■ , (4.5) 

where non-integer powers appear due to the presence of the three -body cut. However, the leading 
non-integer recoil correction at order ;if'^^ dominated by the small scales vanishes, because the 
Pauli principle and the isovector character of the leading nN scattering operator ensure that the 
intermediate A^A^ state is projected onto a P-wave ll40ll . and thus a\ - 0. Specifically, to account 
for the leading correction at order x^^^ in the expression of a"^"' one has to drop small pion 
momenta q with respect to p in the wave functions, which immediately gives zero, as per the 
second line of (14.1b . In consequence the scales '\/M„€, ^M„lm^M„, and yjm^e do not enter at 
this order: all enhanced contributions cancel due to a subtle interplay between the two diagrams 
{d\) and {d^) that is dictated by the Pauli principle. The combined integral is, as originally 
assumed in establishing the ChPT ordering of diagrams, then dominated by momenta of order 
M„. Half-integer corrections at order ;t'^^^ and above can contain the small scales as well, but still 
momenta of order of M„ will have the largest impact on 03 in ( 14.5b . 

The explicit calculation of terms of order x^ and above in this ex pan sion shows that the 
relevant corrections to the static term appear only at orders ;if and x^^'^ 14111 . The net result for 
the recoil correction stemming from these orders is of natural size. In Sect. 14.41 we present the 
results of explicit evaluations of the integrals in (14.1b . 

In principle, there are also contributions with P-wave interactions between nucleons in the 
intermediate state. Examination of the integrand for a,^^ in this case shows that its dominant 
contribution comes from pion momenta |q| ~ ;t'^^lpl' where p is the momentum of the nucleons 
in the intermediate state. Taking |p| ~ M„, we find a contribution from this graph of chiral order 
X^l^ relative to leading. It would appear, then, that this contribution must be calculated explicitly, 
since we need to compute effects at 0(x^^^) in order to achieve our accuracy goal. 

However, this f -wave intermediate-state-interaction graph also includes a factor of the A^A^ 
amplitude, evaluated at an energy e - q^/2M„ — something which was not factored into the 
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Figure 2: Isospin violation in (J3) and (^4). 

above assessment of its chiral order. The A^A^ interaction that appears here must be evaluated 
sub-threshold, but the energy involved is small enough that we can still estimate its effect using, 
e.g., the effective-range expansion. When this is done we find that the A^A^ amplitude will be 
(perturbatively) small at the energies where it is needed for this graph, essentially because the P- 
wave phase shifts are small throughout the domain of validity of the effective-range expansion. 
Hence, we estimate that the contribution to the n^d scattering length of the A^A^ interaction in 
P-waves will be of the order of 



X 



3/2, 



6.p,{E = M2/mp)|fl^'^'-, (4.6) 



where J € {0, 1,2), and, if the phase shifts ^3^^ are small, then the effective-range-expansion 
amplitude has approximately the same size at positive and negative energies of equal magni- 
tude. The A^A^ amplitude thus produces a further reduction beyond the straightforward chiral- 
EFT counting of ~ 0.15, so we find that P-wave intermediate-state A^A^ interactions affect our 
final result by only about 0.2 • lO^^M^', which is significantly below the accuracy that we seek. 

4.2. Further leading-order diagrams 

According to the power counting, the diagrams (c/3) and (d4) (cf. Fig. |2]i should be of the 
same order as double scattering. However, these diagrams are suppressed by accidentally small 
spin-isospin factors 13911 . As the Clebsch-Gordan coefficients are always smaller than 1, it is 
guaranteed that this mechanism can only lead to a suppression and not to an enhancement, such 
that our accuracy estimate remains unaffected. Nevertheless, we have calculated the full isospin- 
violating corrections to this class of diagrams, as required by the power counting, 

..^ SJ^l" (/ q-o-iq-o-2 \ 4Awq-o-iq-o-2 \l 

1287r4^d^^ U (q2 + M^)2 / M^, \ (q^ + M^)^ Ir ^ ' ' 

We find that isospin-breaking corrections due to the pion mass difference amount to 4A„/M\ ss 
28 %. This is a large isospin-violating effect, which is, however, of little practical relevance given 
the overall suppression of this contribution. 



As was shown in 113811 . NLO contributions to the leading-order diagrams vanish in the isospin 
limit. As isospin-breaking corrections to this are suppressed by another two orders, sub-leading 
corrections to idi)-{d4) may therefore be safely ignored. 
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4.3. Triple scattering and the multiple-scattering series 

Although the triple-scattering diagram (see Table [ITl is formally suppressed by p^ relative 
to the double-scattering operator, it was shown in 1381 l39ll that this diagram is enhanced by a 
factor of TT^ over its power-counting estimate and hence has to be included in order to achieve the 
desired accuracy. Neglecting isoscalar contributions as well as isospin-breaking corrections, one 
finds 

TT^d \ Iql / 
The enhancement can be traced back to the occurrence of two successive Coulombic propaga- 
tors. This special topology produces a dimensionless integral that is not 0(1), as assumed in 
dimensional analysis, but instead yields a factor ;r-. This enhancement of the "triangle topol- 
ogy" is not accounted for by Weinberg's application of ChPT counting to the irreducible ttNN 



graphs 113911 . Similar enhancements occur at higher orders in the multiple-scattering series, too, 
and one might worry that this spoils the convergence of the perturbative expansion. However, 
we find from explicit numerical evaluation that the result for the full multiple-scattering series 
resummed in configuration space 156145911 

^ _4 ^ _4(^P« )- r 

differs from the first two terms (double and triple scattering) by only 0.1 ■ 10""' M^'. This is 
significantly below the estimated uncertainty due to the contact term of about 1 ■ lO^^M^' . Hence, 
the multiple-scattering series converges sufficiently quickly that quadruple scattering and higher 
orders can be neglected in the calculation ofa„,j- In particular, we stress that this result indicates 
that the enhancement of triple scattering with respect to its Weinberg-counting estimate does not 
lead to an enhancement of the nNN contact operators contributing to n^d scattering. These still 
occur at O(p^) relative to leading order, preserving the estimates of the contact term given above. 
For a more detailed discussion of the multiple-scattering series in meson-nucleus scattering we 
refer to H. 



4.4. Evaluation of the matrix elements 

Our results for the wave-function averages are shown in Table [3j We give results for two 



different modern phenomenological interactions, AVI 8 116211 and CD-Bonn 1I63I1 . and the im- 
plementations of 1611] of chiral interactions at order N^LO. For the calculations, we used wave 
functions obtained from numerical solutions of the Schrodinger equation in momentum space. In 
order to facilitate the calculation, we used a recently developed Monte Carlo scheme to evaluate 
the integrals. For details on the numerical procedure, we refer to 13^ . The calculation was also 
cross-checked with the standard method of Gaussian numerical integration. 

As the isospin-breaking corrections to the ttN scattering lengths in ( 14.31 ) are relevant only for 
^static^ to which they contribute about 1%, we may write the integrals for the contribution of the 
cut term and the NLO correction to the static ;rA^A^ propagator as 

«™' = %(fl-)2/-', <J5 = ^(fl-)'/N^°, (4.10) 
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model 


A/A 


<l/q-> 


jNLO 


rcut 


<l/lql) 


a"" 


jEM 


AV 18 


— 


12.7 


1.94 


-2.68 


7.4 


-0.00050 


9.92 


CD-Bonn 


— 


12.8 


2.04 


-2.47 


8.7 


-0.00037 


9.85 


ChPT N^LO 


450/500 


13.0 


2.12 


-2.29 


9.7 


-0.00007 


9.85 


ChPT N^LO 


600/500 


12.8 


2.04 


-2.43 


8.8 


-0.00025 


9.94 


ChPT N^LO 


550/600 


12.9 


2.09 


-2.39 


9.5 


-0.00020 


9.91 


ChPT N^LO 


450/700 


13.2 


2.15 


-2.30 


9.9 


-0.00007 


9.95 


ChPT N^LO 


600/700 


12.9 


2.09 


-2.43 


9.4 


-0.00025 


9.85 



Table 3: Matrix elements for phenomenological as well as N^LO chiral wave functions. The cutoffs A/A are given 
in MeV and specify the version of the chiral interaction as given in 16111 . All integrals are given in appropriate powers 
of M^. I^^ is defined in Sect. l5.4l and needed for diagrams involving virtual photons. 



-24.1+0.7 



-•nlo 



3.8 ± 0.2 



-4.8 ± 0.5 



„ triple 



2.6 + 0.5 



-0.2 + 0.3 Afl® 



0.2 



Table 4: Strong contributions to a _ in units of 10 ^^M^ ' for a 



1.1- 10-^ M- 



with 



!^' 



pd\i¥{p)-¥{p-q)mp) 



q2 + 5 q^+6 



2\q^+6 q2 + 5+p/J 



iNLO 



^ / 1 / ^q 

\ q^ ^ Wq + mp 



(4.11) 



and 6, 6, and p defined in ( 14.2b . 

The statistical uncertainties of the numerical results are not significant and are therefore not 
given in the table. An appreciably larger uncertainty is introduced by the different short-distance 
(r <K 1 IMjj) physics of the A^A^ wave functions used. To combine the results for the different 
deuteron wave functions, we take the average of all seven potentials as our mean value, while the 
uncertainty is taken to be the maximum deviation from this average. In this way, we obtain the 
individual contributions to a*^" given in Table |4] (note that Aa*^* is independent of the deuteron 
wave function). They produce a total 



= (-22.6+ 1.1+0.4)- lO^^M^ 



(4.12) 



where the first error refers to the model dependence of the matrix elements and the second to 
the uncertainty in the isospin-breaking shifts in the jiN scattering lengths. (Here, and in Table[3j 
results are quoted for cC - 86.1 ■ 10"^ M^'. When our correlation band is obtained in Fig.|6] 
below, the full cC dependence of a^"^ is taken into account.) 

5. Virtual photons 

The improved Deser formula ( 11.81 ) is derived in an EFT that resums the effects of the photon 
ladder in pionic deuterium. This calculation includes effects due to virtual-photon momenta in 
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the "hadronic-atom regime" |k| ~ r^ with the Bohr radius tb = (aM„) . Our calculation of fl;,-^ 
should therefore include all contributions above this scale. In the ChPT counting where momenta 
are assumed to be of order M^ the leading contributions due to the exchange of (Coulomb) 
photons between the n^ and the proton are given by the diagrams shown in Fig. [3] {d(,), (dj), and 
(cfs)- Photon exchange is perturbative at |k| ~ M„, and the pertinent pieces of these graphs enter 
at O(e^) relative to (di), cf. Tabled Such effects in the other diagrams are of a higher ChPT 
order than considered here. 

However, diagrams (de) and ((ig)-((iio) are reducible in the sense originally defined by Wein- 
berg 13611 . with the ;7rA^A^ intermediate state involving relative momenta ~ yjM^e <s M„. In 
fact, these diagrams are "would-be infrared singular" in the sense that in the limit e — > (and 
for static nucleons) they involve the (singular) matrix elements {Ijc^). This leads to a potential 
enhancement ~ yjM^le for physical values of e. Furthermore, the intermediate A^A^ pair can be 
in an S -wave, such that we must allow for the possibility of A^A^ interactions while the pion is 
"in flight". In the isoscalar partial wave the intermediate A^A^ Green's function will, in particu- 
lar, include the deuteron pole. In that contribution we must separate the low-momentum part of 
these contributions that exactly corresponds to the combined quantum-mechanical effect of the 
Coulomb potential and the n^d scattering length accounted for in the improved Deser formula, 
where the infrared divergence is regulated by the presence of the hadronic-atom binding energy. 

We will now discuss all these effects in some detail. As the appearance of additional scales 
might require modifications of the ChPT counting rules, we will also consider the double- 
scattering diagrams {dg) and ((iio), which formally enter at higher order 

5.7. Diagrams without intermediate-state NN rescattering 

As a first step, we consider the diagrams depicted in Fig. [3] where there is no intermediate- 
state A^A^ interaction. For {d(,) and (d^) we obtain 






a<d,nw ^ r^^ rd^pd^q(fl,-„^"'"(p - q)^(p) + a„-p^'\9yi'{9))- 

(5.1) 

where we have used time-ordered perturbation theory to include the nucleon recoil both in the 
photon and the pion propagator. It is now convenient to split this expression into isovector and 
isoscalar nN interactions. Dropping isospin violation in the scattering lengths, the isospin T = 1 
part becomes 






.(4) ^ ^ d3pd3q(vi''(p) - >I"(p - q))vF(p)-__-____-. (5.2) 

n-^i J |q|(|q| + 6|2oJ^^){^^ + S) 



Similar to the double-scattering process the expansion of {d(,) and (c/g) in M„lm^ induces non- 
analytic terms due to nucleon recoil. However, in contradistinction to a'^"' in ( 14.1b the region 
of small momenta in ( 15.2b is significantly enhanced due to the presence of the photon. As a 
consequence, the expansior|jof ((ie) and {d^) starts from 1/ yjx 



Jdf,)Md!.) ^0 



= — F+^static + ^1 V? + --- • (5-3) 



■^Note that this concerns only the aforementioned expansion in M„lm^, whereas the scale of the individual contribu- 
tions is hidden in the coefficients b,. For example, foo ~ 1/t is dominated by momenta ~ y, which together with 1/ -yf^ 
from j5.3) produces 1 / '^M^e in 15. 6t . 
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Figure 3: Virtual-photon diagrams in n^d scattering. Note that the charges of the particles are shown explicitly to 
illustrate how the photon couples to the charged particles. The effect of mass differences between charged and neutral 
particles is not considered, since it is higher order in isospin breaking. 



Indeed, at leading order in M„/mp, which appears from the non-perturbative regime of the three- 
body cut (Iql ~ -sjM„/mp\p\ <K |p|), the contributions of {d(,) and (c/g) are equal and both involve 
the integral 



f ., .3 *I"(P)*I'(P) 



q2(q2 + 2M„{e + p^/mp)) ^Jx 



I^'"-, 



V2;r2 r., *npmp) 
Vp^ + ^'^ 



with y - -dnipe. If we use asymptotic deuteron wave functions 



*I'(P) = 



Vr 



;r(p2 + y2) 
to perform the integral, we find at leading order 



afj, = -afl\ 



2aM„^pa Stt 1 



(5.4) 



(5.5) 



(5.6) 



In this way, we see that indeed the individual diagrams are enhanced by y/M^/e, but these con- 
tributions exactly cancel in the sum of (c/e) and (c/g), such that ^0 = in ( 15.3b (in close analogy 
to fli =0 in ( I4.5l l). The physical explanation for this cancellation is provided by the Pauli prin- 
ciple: as the ttN interaction does not change the spin, it implies that for the isovector case the 
intermediate A^A^ pair must be in a P-wave, which is reflected by the relative sign between the 



19 



tiO(£iOtt 



Figure 4: Leading virtual-photon diagrams witli tlie black blobs indicating intermediate NN interactions. The time- 
reversed diagrams are not shown explicitly, but are included in the calculation. 

wave functions in ( I5.2l l. As P-wave A^A^ interactions are small (cf. discussion in Sect. 14. Il l, the 
only non-vanishing isovector contribution is therefore generated by the residual effects in ( 15.2b . 
The calculation of the gauged Weinberg-Tomozawa diagram (dj) including the nucleon re- 
coil in the photon propagator proceeds along the same lines as the decomposition of (di) and 
((^2)- Subtracting the appropriate two-body diagram, we obtain: 



len^Fl^Aq^l l67r^Fl^A\q\{2mp + \q\)l 

(5.7) 

where the second and third terms are analogs of fl^LO ^"'^ '^'^"'' respectively. Note that when 
momenta are of order M„ the recoil correction in the photon propagator is, in principle, a higher- 
order effect. Indeed, we find that the corrections to the static photon propagator are numerically 
very small, only about -0.045 ■ lO^^'M^', and may therefore be safely neglected. 

In contrast, the Pauli principle allows for S -wave A^A^ interactions in the isoscalar part of {d(y) 
and (c/g). These will be discussed in Sect. 15.21 while the numerical results are summarized in 
Sect. 15.41 However, if the T = part of these diagrams were significant, then one would also be 
concerned about virtual-photon exchange within the more sizeable, double-scattering, diagram. 
For this reason, we also give the expressions for (c/g) and (c/io). Neglecting the nucleon recoil in 
the photon propagator but keeping it in the relevant ttA^A^ intermediate state, the result reads 



^(rf,)+(rfio) 



it^U \q2 J I2(l2 + 5p,p„)/ it^U \J (q + 1)2 12(12 + Vp+i)i 



(5.8) 



which can be separated into its isoscalar and isovector pieces as follows 
y,)+y,„) _ 2aM^fl-)2 / 1 



^T=Q 



^(dg)+(dm) 



7^*^d 



I I r dh r dH \ 

WJ l2(P+cJp,p„)^J (q + 1)2 12(12 + 5p,p,o/' 

^ J 12(12 + 5p,p„) ~Jiq + 1)2 12(12 + 5p p^,) /■ ^-^ 



Again, the leading, potentially enhanced contribution in the isovector case cancels in accordance 
with the Pauli principle. The isoscalar case, including intermediate-state A^A^ interactions, will 
be addressed in Sect. 
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Figure 5 : Contribution to tlie modified Deser formula in tlie nD atom due to the part of the Coulomb-photon ladder that 
is subtracted in )5.10t and 15. 26) . The double line labels the deuteron, the box con'esponds to n^d scattering, and the 
circle to the photon coupling of the deuteron. 

5.2. The role of rescattering I: single scattering with photon exchange 

The isoscalar contribution of {d(,) and (c/g) including intermediate-state A^A^ interactions (ISI) 
(see Fig. nil is given by 

||i|a=,dV*.(,-,{G.(,'-^,,-|-.-^,-.) 

where Gj(p', p; E, P) is the isoscalar A^A^ Green's function for initial (final) relative momentum 
p (p') and a state of total energy E and momentum P. Note that to obtain this result we have 
neglected the recoil effect in the photon propagator, and treated the pion as a non-relativistic 
particle. In the regime which is the focus of this section, where the pion momentum is much less 
than its mass, both of these are perturbative corrections to the main result. Therefore, we neglect 
them for the time being, and focus our efforts on investigating the infrared enhancement of these 
graphs. The two Green's functions may be interpreted as "direct" and "exchange" contributions, 
i.e. the photon can couple either to the nucleon that undergoes the jiN interaction or to the other 
one. The last term subtracts the deuteron pole {q — > 0^), as this part is already accounted for 
in the modified Deser formula, cf. Fig. |5]for a graphical illustration of this piece. The shift of 
the Is level in pionic deuterium is proportional to the convolution of the n^d scattering length 
with the Coulombic wave function of the atom, which diagrammatically correspond to an infinite 
ladder of Coulomb photons. The simplest example shown in Fig.|5]thus needs to be subtracted 
in (15.10b in order to avoid double counting. The details of the derivation of (15.10b are relegated 
toT 



Appendix B.l 



/ 



The isoscalar propagator Gs is real for energies below the A^A' threshold and constructed out 
of continuum states normalized a^ 

.T^*^p"'"(q)*^P'(q) = (2^)'^*''(P' - P). (5.11) 

with 'Pp(q) obeying 

(P__S_W(q)= r±aly^(q,q')T^(q'), (5.12) 

where Vs is the projection of the A^A^ potential onto the isoscalar part. The free part of G, is 
therefore 

E + ij] - P-^/4m„ - q'/mp 



^Note that this means that our continuum NN wave functions are normalized differently to our deuterium wave 
function. 
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while the total Green's function can be related to the A^A^ scattering T-matrix as 

.(Ov„, „..„^, T{q',q;E,P) 



G,(q',q;£,P) = G*"'(q',q;£,P) + 



{E + ij] - p2/4mp - q2/mp)(£' + ij] - P^/A-nip - q'^/nip) ' 

(5.14) 
where T is connected to the A^A^ phase shifts via 



T{k,k;E,P)^- 



\n 



nip k cot 6(k) — ik 



(5.15) 



with k - Jmp(E — P^/A-nip). Alternatively, Gs can be rewritten in terms of the bound and 
continuum state wave functions 



G,(q',q;£,P) 



E + ir] + e- P^/4-m 



P J (2;r; 



n(q')n'(q) 



)^ E + i7]-p~/mp-P^/4mp' 



(5.16) 



The additional factor (In)^ for the deuteron-pole part is due to our conventions for the deuteron 
wave functions ( I4.4l i. Inserting the free part of the Green's function ( 15.13b into (15.10b reproduces 
the expressions for the structureless diagrams discussed in Sect. 15.11 up to higher-order terms 
neglected in the derivation of (l5.1Qb . Using the decomposition ( 15.16b . the isoscalar contributions 
can be cast into the form 



{df,)+{ds) (df,)+(df,) 

"r=o 7'=o,isi 



where 



l^na^pO' 

j3. 



r d"k \ I |/ 
J (2^k2lk2 



|F(k)|2 - 1 



/2//D - iri 



r TT^e ^ 2/ Lm -G^(k)^(G;'(k) + G;-^(-k)) 



(5.17) 



F(k) = Jd'q^'Hqmq - k/2), G^(k) = Jd3q>P+(q)>P^(q - k/2). (5.18) 

Now, from the normalization condition of *I'(q), and the orthogonality of *P(q) and *I'p(q) for 
vanishing momentum transfer, it follows that 



\F(kf - 1 = O(k'), 



G;(k) = 0(k). 



(5.19) 



In this way, we conclude that also in the isoscalar case no terms enhanced by y/M^/e remain. 
Due to the chiral suppression of a^ the sub-leading corrections can simply be dropped. This 
reasoning used here — based on exploiting orthogonality of bound state and continuum wave 
functions — follows the calculation of recoil corrections to Kd scattering in [41]. 

Alternatively, the cancellation can be derived within HittEFT In this case it is also convenient 
to split the total effect into parts without and parts with an intermediate-state A^A^ interaction. 
The portion without intermediate-state A^A^ interaction gives (cf. ( 15.31 ) and ( 15.6b in Sect. 15. lb 



Jd,) 



Jd6) 



16a 
16^ 




2e \ \ m 

-( 




M„ ^QlMn 



(5.20) 
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Again, these expressions show the anticipated infrared enhancement. 



Using the leading-order form of the A^A^ scattering amplitude in pionless EFT 116446 811 



T,,NN(q',q;E,P) = — , (5.21) 

nip 

y + i^jmp. 



■i^-S;] 



we evaluate the diagrams with ISI and find, after removing the deuteron-pole piece already ac- 
counted for in the modified Deser formula, 

"7-=o,isi - 3;r V 2e' 

which precisely cancels the leading piece of the free part of the diagrams in ( 15.20b . Furthermore, 
any contributions of momenta |k| ~ ^M„e may only appear in the non-integer terms in the 
expansion (15.3b and thus are suppressed by an additional power of M^/nip beyond this (i.e. they 
correspond to the third term in ( 15.3b ). These contributions thus have an overall size 

32a M„ /m; 

— -./——«+ a! 2.8 Q-fl , (5.23) 

37T m-p \ 2e 

and are therefore well beyond the accuracy we claim for our calculation. 

There is still a possible contribution from ISI diagrams with momenta of order y, though. 
This would be enhanced by M„ly compared to its naive ChPT order, and so could be relevant for 
our analysis. Decomposing the Green's function according to 

:(0),„, „. . „^ , (27r)3>P(q')^t(q) 



G,(q',q;£,P) = G^"\q',q;£,P) + 



■ii] + e- p2/4;i 



r7^^(q',q;£,P)- 



£ + ///- PV4OTp - q'2/OTp ■'•'v^'^'^' ' '£ + /77-p2/4mp-q2/OTp' 

(5.24) 

i.e. into its free part, the deuteron pole, and the non-pole isoscalar A^A^ T-matrix 7""^^, 
tional contributions due to intermediate A^A^ interactions are given by 



a 



^="-isi '- J (271)^ kA J (2;r)3 ^"^ '-e- k2/2//o - (q' - k/2)2/m. 



MHd,) ^ _^"€p^^ r d^k If rd^qd^q',^,.^, ,^ 1 

T:% - T7nm{^' - f ' q - I -^ - ^' -k) + T:^n{^' - I -q + |; -. - ^, _k). 



p 

X t:' — ttt; — \^ ■ ,..0 , 't'(q) + 2 '-;;:' ' -. K (5-25) 



.np 1 „„„. .. I^(k)|^-1 



This equation will undergo an explicit numerical evaluation in Sect. 15.41 The part of the integral 
involving momenta |k| ~ M„ only contributes at the naive chiral order of the pertinent diagrams. 
For that effect we only need to evaluate {d(,)-{dg,) without any intermediate-state interaction, and 
it is furthermore sufficient to keep only the T - I pieces of the free parts of these graphs, see ( 15.2b 
and (15.71) above. 
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5.3. The role of rescattering II: double scattering with photon exchange 

If only momenta |k| ~ M„ were important, there would be no need to consider the diagrams 
{dg) and {d\Q). But, the full isoscalar contribution to {dg) and (^lo) reads (for the derivation 
seeT 



Appendix B.2i 



(rf,)+(d,„)^ W,)+Wo) _ 327rV^fl-)2 |- d^k f d^l 1 fd^qdV 



r d"k r An i p 

J '{2n9 J (2^k¥j ■ 



(5.26) 



x>FV){2G.,(q'-|q-^.l;-.-^,-k) 

,2G.(q'-|-q.^-l;-.-^,-k)-li^^)^^2^ 
^2 2 2M;^ -^ -k-/2//o + ;?7 J 

These diagrams, too, are reducible, i.e. would-be infrared divergent. The same arguments as for 
single scattering yield that the leading contributions from momenta |k| ~ y/M^e cancel. This 
leaves an effect from these momenta that has a numerical size a; 1 .4 aajr-d, which is significantly 
below the few percent accuracy that we seek. 

Again, however, we are concerned about momenta of order y, which could yield contribu- 
tions enhanced by M„ly, and so compromise the accuracy of our calculation. To evaluate this 
contribution explicitly note that the free part of the Green's function reproduces (I5.9l l, while the 
intermediate-state A^A^ interactions lead to 



(d9)+(fl'io) 



32n^a{^^a-f r d^^k f d^^l 1 ( fd^^qd^q' ^P'Xq') 



r d'k r dT 1 I r 

J (2^J (2^k¥lJ 



p 



^=0''si ^d J (2;r)3 J (2nf k^P I J {2jif -e - k^/lfio - (q' - k/2)2/m. 

X ..iVA^ _^ _ 1^2/2^^ _ (q _ k/2 + l)2/mp ^ -k2/2;Uo + i^ /' 

f:^. = 2r-,(q'-|q-^.l;-.-^,-k) 

.2ryq'-|-q.^-l;-e-^,-k), (5.27) 

where we have used repeatedly that *l'(p) = ^^(-p). 

5.4. Numerical evaluations 

We explicitly evaluate the isoscalar contributions to single and double scattering for which 
expressions were derived in the previous two subsections. For this purpose we use a separable 
A^A^ interaction, since we anticipate that the integral is dominated by low-momentum modes, so 
details of the potential are not important. Specifically, we use the effective potential 



1 

p2T^2 



V(p,p') ^ Ag(p)g(p'), g(p)^-^—^, (5.28) 



where /I is a constant tuned to reproduce the binding momentum y - -Jnipe, and y6 = 1 .4 fm 'is 



introduced to parameterize the effective range of pn scattering, which enters in realistic potentials 
through the one-pion exchange. Solving the Schrodinger equation ( 15.121 ) for the bound state, one 
obtains the deuteron wave function of the Hulthen type 

>l'(p) = A^4^, N^-Jy/3{y+/3)\ (5.29) 
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Using this wave function and evaluating (15.101 ) with G, from ( 15.141 ). we find 

af^^'-''''^ + af]^^^ = -0.034 a\ (5.30) 

We note that the individual contributions a^^p ' and Ojt^^^^^ are 5-6 times larger than their 
sum, which attests to the cancellation derived in Sect. 15.21 at leading order both Oj'^q and 
'^r=o^isi acquire large contributions from momenta of order '\/M„€ which, however, cancel com- 
pletely in the sum. The deviation from zero in this result is mainly provided by higher-order 
corrections in the expansion of fl;^^Q ^ , according to (15.3b . whereas higher-order corrections to 
Oj-^Qi^i appear to be numerically negligible. 

The dominant effect in this result is ultimately due to momenta of order y, and so should also 
be accessible in heavy-pion effective field theory (HttEFT). From (15.201 ) we see that, in H;7rEFT, 
the free piece of the isoscalar contribution to (de) is 



{de)+{ds) _ _ + ^?r 




-0.022a+, (5.31) 



where the explicitly evaluated contribution arises from momenta of order y. This number is actu- 
ally quite close to ( I5J0I ). which suggests that indeed momenta of order y are largely responsible 
for this contributionjj 

In addition, the ISI part generated by the "form factor" in ( 15.25b can also be calculated ex- 
plicitly within (H;rEFT) 






rd^k 

J 1^ 



4f ' ^' = — T-a— ^m^r - 1). (5.32) 



For momenta |k| ~ y the deuteron wave function can be replaced by its H;rEFT approxima- 
tion (15.5b . which leads to 

Ay \hr\ 

F(k) = -^arctan^. (5.33) 

|k| Ay 

Performing the last integration then yield^l 

^id,)Hd>) ^ _2 ^j ^ jQg4)^^^+ ^ -0.035 fl^ (5.34) 

which is even larger than ( 15.31b . However, the numerical smallness of the full contribution from 
intermediate A^A^ interaction apart from the scale y/M^e suggests that the contribution of the form 
factor should be canceled by the non-pole part of the A^A^ amplitude. 



^Taldng into account the wave-function-renormalization factor Z = 1.690 from 16911 . which is necessary to ensure 
the correct asymptotic 5 -state normalization of the deuteron wave function, 15. 311 changes to -0.037a*, which is even 
closer to )5.30t . 

'strictly spealfing, a* should be replaced by 5* + Aa* in order to account for isospin violation, however this does 
not change the prefactor Moreover, in the expeiimental analysis of the level shift in nD also corrections due to the 
electromagnetic radius of the deuteron are talcen into account. We show in | Appendix C. l| that this amounts only to a tiny 
modification of the prefactor. 
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The numerical evaluation of ( 15.261 1 with the pionful wave functions based on the separable 
A^A^ interaction described above yields 

afX^'^"^ + fl^tojsr' = 0.3 ■ 10-^M, ' . (5.35) 

Again, this result basically stems from residual contributions of the free diagrams from scales 
above y[M^. Contrary to {d^) and (c/g), which are dominated by momenta of order y, in this 
case also momenta of order M„ significantly contribute to the diagrams. On the other hand, 
the "form-factor" contribution Qj-'q ' can still be addressed in HttEFT, along the same lines as 



employed for single scattering (see Appendix C.2 1. The corresponding result for the form-factor 
contribution is large. 



8aM^(fpfl-)^.7r2 . aM„d 
-•PF :r^ ":r~2.yi6 — — 

^ n'-e ^2 '' 2r^d 



^W,)+W,„) _ """^^v^P" > 1'^ _ o Qi ^\ _ '^^^^^" _ n ^^ . 1 n-3 ,,-i 



(^ - 2.916) - ^^i:g— = 0.35 . 10-^M,-', (5.36) 

however, it should again be compensated by the non-pole part of the diagrams, since the residual 
contribution of the A^A^ ISI is numerically negligible. 

These full evaluations show that, despite their ostensible infrared enhancement, the isoscalar 
parts of {df) + (c/g) and (c/g) + (c/io) have contributions from momenta of order y that yield parts 
of Qjj-d which are still significantly smaller than our theoretical uncertainty. We will therefore 
simply drop the isoscalar contributions in what follows. 

Similarly, full evaluation shows that — after the Pauli-principle-enforced cancellation of con- 
tributions of order yfhije — isovector contributions of {dg) + (dw) are very small, only about 
-0.1 ■ lO^^M^'. Thus the same sort of cancellations that preclude the existence of a ;t''^- contri- 
bution due to recoil in the double-scattering diagrams also enforces the smallness of this T = 1 
part of (dg) + (dio). We therefore conclude that (dg) and (dio) may be omitted altogether from 
our analysis. Indeed, had we found that (dg) and (dio) were necessary for a precision evaluation 
of a„-d, we would have been forced to consider all photon diagrams at this order, since (dg) and 
(dio) do not, on their own, form a gauge-invariant set of diagrams. 

In fact, the results found here as regards "would-be infrared-divergent" diagrams are very 
important, as the cancellations we have identified guarantee that the original ChPT power count- 
ing, which assesses the impact of momenta ~ M„ on the integrals, provides a reasonable estimate 
of diagrams involving virtual photons, since the remaining infrared enhancement is too weak to 
severely violate the ChPT estimates. In this way, we are left with the diagrams {d(,)-{ds) 

^ laM^^pQ- _ a I 1 

/™= fd^pd^qCy' (p) - y' (p - q))T(p) ^ . (5.37) 

J lql(lql + 5/2wq)(q-^ + 6) 

With the numerical results from Sect. l4.4l and a" = 86.1 ■ lO^^M,,:', we obtain 

flEM ^ (Q 94 + 0.01) ■ IO^^M;', (5.38) 

where the error again reflects the wave-function dependence as follows from Table [3] Thus, 
virtual photons increase Re a„-d by about 4 %. 

In summary, we have shown that there are no infrared-enhanced photon contributions 
from momenta ~ yfeM^ — due to subtle cancellations both for isoscalar and isovector ttN 
interactions — and that the infrared enhancement provided by momenta ~ y is too weak to gen- 
erate effects that significantly exceed the estimates for momenta ~ M„. The size of the virtual- 
photon corrections is thus roughly in line with the original ChPT power counting. 
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6. Dispersive and A corrections 

There are two additional contributions to the n^d scattering length that have not been men- 
tioned so far. First, diagrams with pure A^A^ or NNy intermediate states yield so-called dispersive 
corrections. It is this class of diagrams that produces the imaginary part of a,j-d, although here 
their leading contribution to the real part of a„-d is suppressed by p^^'^ relative to {d\) 114211 . Di- 
agrams with explicit A degrees of freedom enter at the same order and provide the second class 



of contributions we consider in this section II43I1 . The A(1232) contribution is a true three-body 



effect, since the nucleon recoil is needed if this P-wave resonance is to contribute to S -wave n^d 
scattering. Typical examples for each of these two effects are depicted in Table[T] Both classes 
were computed in [(42, ^^J using a calculation for A^A^ — > dn up to NLO in ChPT 17011 in which 
all integrals were cut off at 1 GeV. We have checked that varying this cutoff does not introduce 
additional uncertainty, and the effect in a„-d is then 

^disp+A ^ (_Q g ^ J 5) . 10-3^;!. (6.1) 

Since this contribution is only 0{p^^^) larger than the contact term, we need not include isospin- 
violating corrections to 0'*'^?+^, which, counting e ~ p, would be suppressed by another two 
orders. 



7. Summary of three-body contributions to the pion-deuteron scattering length 

We now summarize the content of Sects. |3}{6]by listing the three-body corrections that need 
to be taken into account in the actual calculation. We recapitulate the numerical results for these 
contributions, and analyze the different sources of uncertainty in the calculation. 

The essence of the discussion in Sect. |3] was to demonstrate that all isospin-conserving three- 
body corrections can be reliably calculated up to 0(p^^^), i.e. a relative accuracy of (M^/nip)^^^. 
This is half an order lower than the contribution of the leading unknown {N^N)^nn contact term, 
which is 0{p^). The uncertainty anticipated due to the truncation of higher-order terms is a few 
percent, as follows from a naive dimensional analysis. However, to achieve this accuracy one also 
has to account for three-body isospin-violating corrections, which are suppressed by e^ compared 
to the leading isospin-conserving terms. Therefore we also presented a complete calculation of 
the isospin-violating corrections up to 0(e^p^^^). To this order the three-body isospin-violating 
corrections are of electromagnetic nature, that is they appear due to virtual photons and the pion 
mass difference. The quark mass difference starts to contribute only at 0{e^p). 

In order to reach our accuracy goal, in Sect. |4] we explicitly evaluated the diagrams of a 
multiple-scattering topology as well as those that involve 3nN^'N and An vertices (cf. Fig.[T]). In 
particular, in Sect. 14. II we derived the expressions for the most relevant double-scattering opera- 
tor including isospin-violating corrections. To better understand the relevant scales contributing 
to the process the full double-scattering expression was divided, without making any approxi- 
mations, into three parts: the static term at LO, its correction at NLO, and the contribution that 
contains the three-body cut. As stated above, the leading isospin-violating correction is due 
to the pion mass difference and thus it appears explicitly in the three-body-cut contribution with 
charge exchange, see ( 14.11 ). In order to have the n'^nn and ir^pn thresholds at the proper places we 
also kept the nucleon mass difference in the propagators, although this is formally a higher-order 
effect {0{e~p)). Due to the presence of the cut the effect of the nucleon mass difference could 
be enhanced. In fact, it proves to be negligibly small: the direct evaluation gives just a 0.5 % 
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correction to the static term compared, e.g., to a 2 % correction from the pion mass difference. 
In this way, we conclude that the effect of the quark mass difference in three-body operators is 
negUgible. 

Other diagrams of a multiple-scattering topology appear at 0{p) relative to the double- 
scattering diagram, according to Table [T] The effect of the nucleon recoil in these diagrams 
can thus be neglected there, with, e.g., cancellations enforced by the Pauli principle guaranteeing 
that the nucleon recoil in the triple-scattering term starts to contribute only at 0{p^). Thus, it 
suffices to evaluate the triple-scattering diagram in the static approximation, see ( 14.81 ). The effect 
of all multiple-scattering diagrams beyond triple-scattering was also evaluated in the static limit 
and shown to be negligible, cf. Sect. 14.31 Similarly, isospin-violating corrections are irrelevant 
for all terms in the multiple-scattering series beyond double scattering. 

The amplitudes {d-i) and {dn), which involve 3nN^'N and An vertices, were calculated in 
Sect. 14. 2l including the leading isospin-violating corrections due to the pion mass difference. The 
contribution of these diagrams appears to be suppressed numerically due to accidentally small 
spin-isospin factors 13911 . 

The combination of all effects discussed thus far in this section defines the "strong" contri- 
bution to the n^d scattering length a*^" . The numerical result for a^" was presented in Sect. 14.41 
via evaluations of the wave-function averages using modern phenomenological A^A^ potentials, 
AVI 8 fill and CD-Bonn fHl], as well as chiral A^A^ interactions at order N^LO |6l!l. We found 
(withfl- = 86.1 • 10-^M^i) 

a'"" = (-22.6+ 1.1+0.4)- 10"^ M;', (7.1) 

where the first uncertainty arises from the different short-distance (r «: 1 /M„) physics of the 
A^A^ wave functions, and the second from the isospin-breaking shifts in the ttN scattering lengths. 
The variation in the results due to the use of different wave functions is about 5 %. This provides 
an independent confirmation of the contact term's effect. 

The combined effect of the dispersive corrections and the A(1232) contributions at 0(p^^^) 



was discussed in Sect.|6]and taken from 1421 143 1 



^disp+A ^ (_Q_g +15). iQ-^M^ ' . (7.2) 

Isospin-violating corrections to these diagrams are irrelevant at the order to which we work. 

Finally, Sect. |5]was devoted to a thorough investigation of the effects related to virtual pho- 
tons. Due to the presence of photon and pion propagators these diagrams are infrared enhanced. 
Therefore, keeping the full dynamical structure of the ttA^A^ propagator (including the nucleon 
recoil) is mandatory to avoid infrared-divergent integrals. To the order we are working, it is suf- 
ficient to consider the diagrams {d(,)-{di<,), which form a gauge-invariant set of diagrams at order 
O(e^). The explicit expressions for these diagrams were derived using time-ordered perturbation 
theory, cf. ( 15.37b . Note that for the gauged Weinberg-Tomozawa diagram {d^) no ttNN propaga- 
tors emerge. As a consequence, that contribution is infrared finite even in the static limit, cf. the 
second term in the first line of ( 15.37b . The inclusion of the recoil in the photon propagator was 
also considered, cf. ( 15.7b . It is, however, a negligible effect. The ultimate result for diagrams 
{d(,)-(ds) is then 

flEM ^ (0.94 + 0.01) ■ lO^^^M^', (7.3) 

with only a ~ 1 % wave-function dependence (see last column of Table O. Motivated by the 
large magnitude of the double-scattering term, we also investigated associated virtual-photon 
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Figure 6: Combined constraints in the a* -a plane from data on the width and energy shift of nli, as well as the nD 
energy shift. Figure from (l|]. 



corrections, cf. diagrams (c/g) and ((iio), although formally they contribute only at 0{e^p^) in the 
power counting. The explicit computation of these diagrams showed that the magnitude of these 
corrections is far beyond the accuracy we seek. 

The three pieces a^", a'^'^p^'^, and a^^, when added together, constitute the three-body con- 
tribution to the n^d scattering length. In fact, to a large extent, the novel three-body effects 
computed in this study accidentally cancel 



AflP) 



+ a 



NLO 



. fl™' + flEM ^ (Q_^ ^ Q 7) . 10-3^^1. 



(7.4) 



This cancellation is, in itself, somewhat remarkable, since, e.g. a^LO ^^ ~ ^^ times larger than 
the final central value. The effect of the cancellation is that the main impact of our analysis on 
the extraction of pion-nucleon scattering lengths is our consideration of NLO isospin-breaking 
corrections — in particular the large shift Aa^ = (-3.3±0.3)-10"^Mj^' — in the ttA^ amplitude 17111 . 



8. Pion-nucleon scattering lengtlis 

Combining the dependence of the n^d scattering length on 5^ and a" and the results for jiH 
discussed above, we find the constraints depicted in Fig.|6] The combined lo" error ellipse yields 



(1.9±0.8)-10"^^M^', 



(86.1 ±0.9)- lO-^M-\ 



(8.1) 



with a correlation coefficient pa-a+ - -0.21. We find that the inclusion of the nD energy shift 
reduces the uncertainty of 5^ by more than a factor of 2 and the correlation between 5^ and a' 
by more than a factor of 3. Note that in the case of the nH level shift the width of the band is 



29 



eP Aa ,Aa^^-^ ^3'^^^^"^'^ Aa^ Wave-function averages 
16% 21% 75% 30% 53% 

Table 5: Individual contributions to the en'or on 5^ are added in quadrature to obtain the uncertainty depicted in the 
bands of Fig. [6] The impact of each source of error is given as a percentage of the total (where the second column gives 
the uncertainty in the isospin-breaking shifts ofnN scattering lengths that occur in a"^'' , cf. J4.3t ). Table from Q|]. 

dominated by the theoretical uncertainty in tia„-p, whereas for the nH width the experimental 
error is about 50 % larger than the theoretical one. 

Table |5] shows the individual contributions to the nD error band: as with the tiH level shift, 
the experimental error is much smaller than the combined theoretical uncertainty, whose largest 
individual contribution is produced by the uncertainty in a^isp+A -pj^g wave-function averages 
contribute about 0.5 ■ lO^^M^ ' to the overall uncertainty in 5+, which is in line with the estimated 
impact on a,j-4 of the 0{p^) — relative to (c/i) — contact term. 

To deduce a value for a^ itself, further input on c\ and f\ is required according to (12.4b . 
c\ is related to the nN cr-term: (r„N - (45 ± 8)MeV as quoted in 111 111 corresponds to c\ — 
(-0.9 ± 0.1) GeV"' and ci = (-1.0 + 0.2) GeV"' at third and forth chiral order, respectively ||3- 
Recent determinations based on nN threshold parameters yield c\ - (-0.93 ± 0.07) GeV ' 12611 
andci — (-1.2 + 0.3)GeV"' il5ll . while an investigation oinN scattering inside the Mandelstam 
triangle led to ci = (-0.81 + 0.12) GeV ' III- Finally, fits of ChPT amplitudes to phase-shift 
analyses provide values in the range ci - -(1.2 . . . 1.4)GeV"^ jlOll . In conclusion, we consider 

ci =(-1.0 ±0.3) GeV"' (8.2) 

as a reasonable average of the present knowledge on this LEG. Taken together with the rough 
estimate |/i| < 1.4 GeV"' 176117711 . this value for ci and ( 18.11 ) yield a non-zero a^ at better than 
the 95 % confidence level 

fl+ = (7.6 + 3.1)-10"^M^'. (8.3) 

The final result for a^ is only a little larger than several of the contributions considered in this 
work. This emphasizes the importance of a systematic ordering scheme, and a careful treatment 
of isospin violation and three-body dynamics. A reduction of the theoretical uncertainty beyond 
that of the present analysis will be hard to achieve without additional QCD input that helps pin 
down the unknown contact-term contributions in both the nN and nNN sectors. 

Finally, we can combine our values for the scattering lengths in the isospin limit with the 



isospin-breaking corrections 112411 to arrive at the nN scattering lengths for the physical channels 
summarized in Table |6] Note that the difference between scattering lengths in the same isospin 
channel is better known than the scattering lengths individually, since the scattering lengths in 
the isospin limit and the associated uncertainties drop out. For example, the difference between 
a„Op and a„o„ at NLO is given by 12411 



^Tfip - CIrfin - 



1 1 4c5B(md - mu) Ml , , , , 

W} f] " -^(VA. + 2M.A.-VA.-2M.A.) 

(-3.4 + 0.4)- 10"^^ M^', (8.4) 



'^This value is consistent with recent determinations of cTj^n from the lattice, (T^n = (50 ± 10)MeV 1721 and cr,^ 
(50 ± 10 ± 10)MeV01. 
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isospin limit 


channel 


scattering length 


channel 


scattering length 


a+ +fl" 


n^ p — > n^ p 


86.1 ± 1.8 


n^n- 


-> n^n 


85.2 ±1.8 


a*-a- 


n^ p — » TT^p 


-88.1 + 1.8 


TT^n - 


-> n^n 


-89.0 ±1.8 


-V2a- 


TT^p -^ jfin 


-121.4±1.6 


K^n- 


■^n^p 


-119.5 ±1.6 


fl+ 


n'^p ^ ttV 


2.1+3.1 


n%- 


^n% 


5.5 ±3.1 



Table 6: nN scattering lengths for the physical channels in units of 10 ^M^^ ' . 

with C5 being related to the strong contribution to the proton-neutron mass difference. Although 
it is formally of higher chiral order, here the contribution from the cusp due to n^n and n^p 
intermediate states has been kept, since it is enhanced by half an order in the isospin-breaking 
parameter 6 - {e^,md - niu). And indeed, it ultimately contributes about 30% to the number 
quoted in ( 18.41 ). 

9. Goldberger-Miyazawa-Oehme sum rule 



The Goldberger-Miyazawa-Oehme (GMO) sum rule 111211 relates the charged-pion-nucleon 
coupling constant g^/4-n to a^-,, - a„+p and the integral over the cross sections crj^- „(^) - <^'^+pil<^) 
measured in the laboratory frame for pion momenta k ranging from zero to infinity 

(mp+mn\2 u/ M„\M„ Ml 



l^((^)'-')((-a>— >-t4 



1 r .,<°V(^)-<°V(^) 



r^-r^\ dk '"' ■" —. (9.1) 

This result is derived by writing down dispersion relations for n^p — > n^p for fixed f, assuming 
that the amplitudes are analytic functions of s with a right-hand cut starting at ithr — (m^ + M„)^ 
and a left-hand cut starting at i = (mp - M„)- - t. These dispersion relations are then evaluated 
at threshold. The scattering lengths enter as the ttN amplitude at threshold, while the coupling 
constant gc is related to the residue of the nucleon pole. Finally, one employs the optical theorem 
in the laboratory frame to replace the imaginary part of the amplitude by the total cross section. 

9.1. Isospin violation 

There are two ways in which isospin violation affects the derivation of the GMO sum rule: 
mass effects and virtual photons. 

The proton-neutron mass difference enters through the intermediate neutron in the nucleon- 
pole diagram, which is akeady taken into account in ( 19.11 ). Additionally, the threshold for n^p — > 
7r"« lies below Sthr, such that the right-hand cut for n^p already starts at {m^ + M„o)^. Thus, the 
total cross section for n^p diverges at threshold due to the lower threshold of n'^n. However, this 
divergence corresponds just to the right half of a principal-value integral: the dispersion integral 
for the reaction n^p really starts at i = (»!„ + M^)^, and the resulting pole at Sthr can be taken 
care of in the usual way by the principal-value prescription. To estimate the remaining effect, 
one may use the fact that the imaginary part of the amplitude for s < (wtp + M„-)~ can be well 
approximated by the imaginary part of the n^p scattering length due to the n^^n intermediate 
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state 112411 . In this way, we find a shift in J by about -0.005 mb, which we will take into account 
in our uncertainty estimate for J^ below. In conclusion, mass effects do not invalidate the GMO 
sum rule and the necessary modifications are quite well under control. 

In contrast, we cannot write down the GMO sum rule in the presence of virtual photons, as 
e.g. the nucleon pole is not separated any more from the jN cut. Therefore, ( 19.11 ) is only applica- 
ble if all ingredients are purified from virtual-photon effects to ensure that the analytic structure 
of the corresponding amplitude coincides with what was assumed in the original derivation. For 
this reason, we will adopt the following point of view; we assume that the removal of electromag- 



netic effects in the cross sections using the Tromborg procedure II78I1 works sufficiently well that 
the resulting value for J is compatible with the above analyticity assumptions. Moreover, we 
subtract virtual-photon effects in the scattering lengths based on I24j1 . but keep isospin-violating 
effects due to the nucleon and pion mass difference (we will dwell on this procedure in Sect. 19.31 ). 
In this way, our final result for gc consistently refers to the scenario where all particle masses are 
fixed at their physical value, but virtual photons are switched off. 

Finally, we comment on the definition of scattering lengths and coupling constants in the 
presence of electromagnetic interactions. Even in principle, the calculation of electromagnetic 
corrections is a scale-dependent procedure 117911 . which, however, can only be systematically 
addressed if the underlying theory is known. Within an effective theory, ChPT in our case, a 
consistent treatment of electromagnetic corrections is possible, apart from the fact that the am- 
biguities in the separation of photon effects present in full QCD should be reflected in additional 
uncertainties in the LECs. To the best of our knowledge, the practical consequences of 17911 for 



ChPT calculations have yet to be explored. However, the study of the linear cr model in 117911 
suggests that such effects are not relevant at the level of accuracy at which the LECs can be 
usually pinned down. In addition, the definition of a scattering length for charged particles is a 



subtle matter 118011 . and recent attempts to define a strong proton-proton scattering length yield 



tempts 

Hl-ii 



only scale-dependent quantities Il81l - l84tl . An effect analogous to that discussed in these works is 
also present in the calculation of the nN scattering lengths, but, due to the perturbative nature of 
nN dynamics, this effect is negligible as we will show in the following. Removing the Coulomb 
phase 0c(lpl); the behavior of the n^p scattering amplitude T„-p at threshold is given by 11711 

^-2,-.«c(lpl) T^_^ ^ ^ + B2 log ^ + r*; + O(lpl), (9.2) 

IPI tJ-H 

where 

An^^a„-p = r*'^, Bi = A7i'^aM„a,-p, B2 = -SnaM^ia^pf. (9.3) 

The scale ambiguity represented by the presence of log/i/f is not induced by an ultraviolet di- 
vergence, but by a kinematic singularity at threshold. While the l/|p| term (the leading approx- 
imation to the Gamow-Sommerfeld factor Ii85u86il ) can be unambiguously separated, the log |p| 
requires the choice of a scale in order to define the strong threshold amplitude T^^ . In ( 19.21 ) that 
scale has been chosen to be jU//. However, B2 differs from zero only at two-loop level, i.e. it is 
suppressed by two chiral orders compared to the accuracy at which the isospin-breaking correc- 



tions ll24tl are known. Thus, choosing the mass of the p meson, rather than the reduced mass fiH, 



shifts a„-p by 



-^ log ^ = -2a^iH{a,-py log -^ = -0.2 ■ lO-^M,', (9.4) 



3P 

an effect fully in line with its two-loop estimate that can therefore safely be neglected 



32 



One might worry that our definition of gc is not exactly what is measured in experiment, be- 
cause, despite the appHcation of electromagnetic corrections, the full range of virtual-photon ef- 
fects is not captured by present-day analyses. For this reason, one could try to add a certain class 
of virtual-photon diagrams in order to obtain a quantity that corresponds better to the experimen- 
tally accessible one. In fact, this is quite a difficult enterprise: to extract the coupling constant, 
we need the amplitude at s = m^, where threshold ambiguities as in the case of the scattering 
lengths do not occur However, the cancellation of infrared divergences that is ensured at thresh- 
old by phase-space arguments no longer works, which makes the inclusion of bremsstrahlung 
inevitable. In order to estimate the size of such effects, one may in a first rough approximation 
consider only the leading bremsstrahlung contribution that involves logarithms of the detector 



resolution Zsmax- In this naive approach — described in more detail in Appendix D.2 — we find a 
shift of about 0.07 in g^/A-n for iSmax - lOMeV, which is thus significantly below the accuracy 
we claim for our final result below. We conclude that, in order to address virtual-photon effects 
systematically, one is forced to perform the full radiative corrections for a given process, which 
is beyond the scope of this work. 

9.2. Evaluation of J^ 



The evaluation of J has recently been discussed in great detail in 111311 and 111411 . hereafter 
referred to as ELT and AMS. The main difference between both analyses is that the former 
relies on phase-shift solutions to determine the cross sections, while the latter uses data directly. 
Both investigations apply the Tromborg procedure I78i1 to remove electromagnetic effects|3 The 
quoted results 

j:^^^ = (-1.083 + 0.032) mb, /^ms = (-1-060 ± 0.030) mb, (9.5) 

are consistent within the uncertainties and we take the discrepancy between the mean values as 
an indication of the final accuracy one can hope to achieve in an evaluation of J . In order to 
obtain an average value of J that combines ELT and AMS, we now compare both analyses in 
the low-momentum {k < 2GeV/c), the high-momentum (2GeV/c < k < 240GeV/c), and the 
Regge regime [k > 240 GeV/c) separately. In general, we employ the uncertainties quoted by 
ELT, whose error estimates tend to be more conservative than those of AMS. 



In the low-momentum region, we average the mean of the results for the SM95 1I13L 188 1 



SM99 li3l|89|], and FA02 ||90l|91I] phase-shift solutions with AMS. In these determinations, the 
threshold region constitutes an additional source of uncertainty due to a lack of very low-energy 
data. Therefore, for the interval (0-80) MeV/c an interpolation between the cross sections at 
threshold 

Acr(O) = cr^°'^,(0) - <°',/0) = 4n{{a„-„f + (a^p^ _ (^^^^^^2^ (9.5) 

and the available scattering data is needed. To estimate the impact of the scattering length on J, 
we consider the S -wave part of this interpolation, which changes by about 0.009 mb if one varies 
Acr(O) by 30 %. In view of the fact that the difference between our value for Acr(O) and those 
of AMS and ELT lies below 20 %, this should provide a conservative estimate of the additional 
uncertainties to be expected in the threshold region. Adding in quadrature this estimation, the 
uncertainty quoted by ELT, and the effect due to the lower Tr'^n threshold discussed in Sect. 19.11 



'Above the energy range where the Tromborg corrections are available the effect due to the Coulomb barrier is 
accounted for following the potential-model calculation 18711 . 
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k in GeV/c 


0-2 


SM95 [13, 88J 


SM99 |13, 89J 




FA02 [90, 91J 




-1.302 


-1.314 




-1.3043 




AMS 


average 








-1.3003 


-1.304 ±0.023 






2-240 


ELT[13j 


AMS [14] 




average 




0.197 


0.2149 




0.206 ± 0.024 


240-00 


Hohler [92j 


Donnachie-Landshoff 


m} 


Gauron-Nicolescu [941 




0.0222 


0.0294 




0.0244 




Regge94 [95] 


Regge98 [96j 




average 




0.030 


0.018 




0.025 ± 0.007 



Table 7: Contributions to / from different momentum regions in mb. For tlie detailed comparison, we use that the 
regions (2-2.03) GeV/c and (240-350) GeV/c yield a contribution of 0.0027 mb and 0.0043 mb to r , respectively Q. 



yields the +0.023 error given in Table|7] In the high-momentum region we use the mean of AMS 
and ELT with the ELT error, while the contributions from the Regge regime above 240 GeV/c 
are determined as the average of the five models employed in AMS and ELT with an error chosen 
generously to encompass all models (cf. Table|7]i. In this way, we finally obtain for our average 



J =(-1.073 + 0.034) mb. 



(9.7) 



9.3. Results for the pion-nucleon coupling constant 



Inspired by |13|], we adopt the following strategy to determine a„-p - a„+p. Writing 



^n p 



■ a„^p = 2a„-p - {a„-p + a„-„) + a„-„ - a„*p = 2a„-p - 2(5+ + Aa+) + a„-„ - a„^p, (9.8) 



we can take 



-^n p 



(85.66 + 0.14)- IO^^^M;' 



directly from the level shift in ttH, 5+ from (18.1b . and 



(9.9) 



Qj^ /7 ^n^p 



An^, 



[f2+2M„{2gl+g\)] 



(9.10) 



from 12411 . As the g\ cancel between a„-„ - a„*p and 2Afl+ (cf. (I2.6l l and (19.10b ). such that only 
the rather well-determined LEC /2 remains, this procedure is particularly stable with respect to 
unknown LECs. We find 



3a^-1 



fl.-p-fl.> = (173.2+ 1.6) -IO-^M; 



(9.11) 



However, these scattering lengths still contain virtual-photon effects, such that we need to sub- 
tract the corresponding contribution (cf. [Appendix D.l"| 



-^K p aj^+p\ 



(2.1±1.8)-10"^M^'. 



(9.12) 
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The GMO sum rule ( l9Tt with the input ( |977] |. ( |9^ . and ( |9^ then yields 

;?^ 

— = 13.69 + 0.12 + 0.15, (9.13) 

An 

where the first error gives the uncertainty due to the scattering lengths and the second that due 
to J. This value is in agreement with determinations from A^A^ {gl/A-n — 13.54 + 0.05 19711 ') 
and nN {gljAn = 13.75 + 0.10 IH, gljAn = 13.76 + 0.01 Q) scattering data. We stress that 



the errors quoted in i9lll97ll98ll mainly reflect statistical uncertainties. The systematic subtleties 
associated with isospin violation that were discussed above in Sect. 19.11 were not quantified in 
these previous studies! '°l 

10. Conclusions 

We have presented a calculation of the n^d scattering length based on ChPT including full 
isospin-violating corrections in both the two- and three-body sectors. In the isospin-conserving 
part we included all contributions below the order at which an unknown {N^N)^7r7r counter term 
enters. The remaining cutoff dependence of our results is consistent with dimensional analysis of 
the counter term, and thus provides an estimate of the theoretical uncertainty of the calculation. 
We have also considered isospin-violating three-body contributions below the order 0{e^p) rel- 
ative to the leading isospin-conserving operator Moreover, we considered several higher-order 
diagrams that could potentially be enhanced due to the dominance of double scattering in a„-d 
or because of small scales associated with the deuteron binding energy. In the end, we find no 
significantly enhanced virtual-photon effects. Ostensibly enhanced contributions from momenta 
of order '\/M„e vanish for both isovector and isoscalar ttN scattering, where the cancellation can 
be traced back to the PauU principle and the orthogonality of deuteron and continuum wave func- 
tions, respectively. We conclude that a higher accuracy in both the two-body isospin-violating 
corrections and the three-body part of the n^d scattering length requires additional information 
on nN and nNN contact terms, which are not constrained by chiral symmetry. Finally, we pre- 
sented the consequences of the calculation for the nN scattering lengths and — after carefully 
revisiting the GMO sum rule in the presence of isospin violation — for the charged-pion-nucleon 
coupling constant. 
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Appendix A. Effective Lagrangians 

For the sake of completeness, we review here the effective Lagrangian for nucleons, pions, 
and virtual photons, as constructed in 12611 . The following terms are needed in the context of the 
present study 

£f' + 4"' + Xy = ^{d'^uU.u +/U + U'x) + zf\ququ^) - \f,,f''' - ^{d.A^'f, 
+ k3{{d^'U^QU}{d^.U^QU} + {d^UQU''){d^UQU''))\, 

4^'' = ^^'i'{d,ix-,u,]iy}^ + h.c., 
4'' = F^^fiiQl - 2-) +/2<e+)e+}*i', 

+ £^*i'{^6<e+><e-«,,>Z>' + gi{Q^u,)Q^D>' + g,{Q-u,}Q^D^')^ + h.C, (A.l) 

where (A) denotes the trace of a matrix A, A - A - {A}/2 its traceless part, '^{O + h.c.)*}* = 
'VO^ + h.C. for an operator O and 

di,U ^d^U-iAi,[Q,U], ;tr = 2Bdiag(OTu,md), t/ = M^ S = ^ diag(2,-l), 

F^y ^ d^Ay - dyAf,, g = g diag( 1 , 0) , Q^ ^ -(uQu^' + u'^ Qu) , 

D^^dft+ F^, r^ = ::{uHd^ - iQA^)u + u(d^, - /gA^)M"'"), x± = "'-^m' ± ux^u, 
Uf, = i(u'(dfj - iQAfi)u - u{d^i - /2A^)m' j, [D^,Mv] = df^Uy + [r^,Mv]. (A. 2) 

*P = {p, nY contains the nucleon fields and the matrix U collects the pion fields. F, g, and m 
are the pion decay constant, the axial charge, and the mass of the nucleon in the chiral limit, 
respectively. The renormalized LECs are denoted by a superscript r. 
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Appendix B. Photon diagrams in chiral effective theory 

Appendix B.l. Single scattering with photon exchange 

To begin with we discuss the case that the nNN -^ nNN operator contains only the (nom- 
inally) leading isoscalar contribution. Writing the (isoscalar) threshold nN amplitude as T^ - 
An^pa'^, we have 

J (271)4 J (2;r)3 ^^' (M„v + ky- - Ml + ir] "'^I'k^ 



[iGs{q' 



k k \ / k k \^ 

y>q - y; -e - ^'0, -k) + iG,[q' - -,-q+ -; -e - ko, -k)k'(q), 



(B.l) 



where v = (1,0) and we have already used the fact that ^o ~ ^'/M„ <K k and therefore can be 
neglected in the photon propagator. The factor of 2 multiplying T^ is present because the pion 
can interact with either the neutron or the proton, while the overall factor of 2 includes the time- 
reversed diagram. The factor e/2 in the first line occurs because we include only one of the two 
possible interactions of the nucleons inside the nucleus with the photon. The other is accounted 
for via the exchange term, which is represented by the second Green's function. In that portion 
we have replaced p by -p as compared to the direct piece of the amplitude. More specifically, 
we can rewrite Ai as 

M = (>I'ieG(£)i(l - PvJT^nI'V), (B.2) 

where 

T„-N = (r+1*" + T-rf) (8) 1*2' +(1^2), e = - (l"> + T^,") (8> 1*^' + (1 ^ 2), (B.3) 

the superscript referring to nucleon 1 and 2, respectively. In ( IB.2b G{E) = 1 /{E + irj - H) denotes 
the Green's function describing the propagation of the A^A^ pair from the nNN interaction to the 
photon coupling, and the projector (1 - Pi2)/2 has been introduced to impose the Pauli principle 
(f 12 interchanges nucleons 1 and 2). Note that the operators in round brackets in (IB. 3b are 
implicitly understood to be accompanied by "shift operators" S^'^ indicating the momentum shift 
induced by the pion-nucleon or photon-nucleon interaction, which has to be taken into account 
when the symmetry properties of the individual terms are analyzed. Inserting (IB. 3b into ( IB.2I ). 
we obtain two distinct contributions (all terms with a single Tj involve {T - OIt^ |r = 0) = 
and may thus be dropped). First, the isoscalar piece reads 

m^G(l - Pn)2T^m = {"¥1^(0, + G,,)2T^m, (B.4) 

and since the isospin wave function of the isospin-zero state is akeady antisymmetric under 
particle exchange, this isoscalar "direct-minus-exchange" G(l - P12) = Go - Ge prescription 
(in ( IB .41 ) Gs denotes the outcome for the exchange part) produces the sum of two G,, terms 
in ( IB. lb . Second, the "direct-minus-exchange" contribution to A1 that results from the isovec- 
tor part of the pion's interaction with the A^A^ system discussed below only contributes for 
odd partial waves due to the Pauli principle. This can already be seen from the isospin struc- 
ture: TJ'\T - 0} oc\T - 1,T3 - 0), and since spin is conserved, the Pauli principle in the form 
|-_ j-jL+s+r _ _Y requires an odd partial wave. 
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Because of the presence of the deuteron pole, the expression dB.lb is infrared divergent (not 
just "would-be infrared divergent")- It blows up because we have neglected effects due to the 
atomic-binding energy B^^: if we included these, it would be regulated at scale V^WJAt ~ aM„. 
However, the physics associated with this momentum scale was already included in the atomic- 
physics calculation. In that calculation, which is done using non-relativistic effective field theory 
(NREFT), we have an expression for these effects that corresponds to a structureless deuteron. 
Using the fact that the normalization of any deuteron wave function equals 1, that expression can 
be written as 



, r d'^k r 

J (2^J 



'^^iR=2 — — d\d\'¥(q')i2T^nq') 



-k() - k2/2md + IT] ko - \i?-l2M„ + irj 



X ^\(i){-ie)(ie)^^'(q). (B.5) 

From the difference of dB.ll ) and ( IB.5I ), we can obtain an expression that is safe in the infrared 
and includes only the effects not already accounted for in the NREFT computation 

^iRsafe - 2^^ j (2;r)3k2J {2nf ^^"^Aj 27r T '^'^ 2'** 2' ^ "'' ^1 



k k 



iGs{t( - y, -q + t; -e - ko, -k) 



2M„ 2(2n)Mq')'i'Hq) 



(M^v + fe)2 - M2 + ir] -k2/2;UD + ii] 

(B.6) 



Now, since we are prepared to ignore ttttNN cuts (they only lead to higher-order effects), we 
drop the "backward-going pion" contribution. The evaluation of the ^o integral in the first term 
can then be done by picking up the pion pole. This yields 



where we have neglected terms that are higher order in k, and so replaced Wk by Mjr and the 
kinetic energy of the pion by its non-relativistic form. From this result, we can now read off ( IS.lOl l 
via 

(d(.)+{di) _ ^^hR safe ,r, ox 

Finally, the isovector P-wave part produces a contribution corresponding to 

3i 



_8^ rdk r,3 ,3,>pt(q, 
=' {2nf^i J k2 J ^ ^ '^ 



\oM-\ 



q _ !i; _f _ J^, _k) - gU - -, -q + -; -f - — , -k)W(q), 
^ 2 2M„ ' ^^ 2 ^ 2 2M„ '] ^ 

(B.9) 



where the isovector Green's function G,, is defined analogously to G, with the deuteron pole 
removed and r"^^ — > T"^^ (the relative sign is due to the symmetry of the \T - 1 , 7^3 - 0) 
state). Although we do not consider isovector A^A^ interactions, it is a valuable check of the 
calculation that the free part of the Green's function reproduces the expressions of Sect. 15. II 
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Appendix B.2. Double scattering with photon exchange 

We now generalize the previous discussion to a nNN kernel given by double jiN scattering. 
Let us first consider the case where the photon is exchanged after the nNN interaction. Two 
distinct processes contribute to the nNN interaction: the exchange of a tt" corresponds to double 
elastic n^N scattering, while an intermediate n^ requires two charge-exchange reactions. The 
amplitudes can be written as 

(B.IO) 
where 

T„-M^^N = (- V2 r-)T*_'> ® 1'2) + (1 ^ 2), 

Tn-'N^.-N = (- V2 r-)T*_'> ® 1<2) + (1 « 2), (B.ll) 

with raising and lowering operators 



^ 2 



^(rf±/r^'), (B.12) 



and T„-N and Q are defined in ( IB.3b . Again, the momentum shifts are understood to be taken 
into account implicitly. The evaluation of A^'^ proceeds in close analogy to (IB. 2b . once the con- 
tributions for which both nN interactions happen to the same nucleon are excluded. Neglecting 
isoscalar nN scattering terms (which are of higher order), we obtain 

M"' = -2(T-fm^G^il - Pn)m = -2(r-)2(v}'|£(G, + G.W). (B.13) 

The expression for Al'^ is somewhat more involved, the crucial observation being that the states 
^(2)^U)|7. = 0), T^^W^^^IT = 0), (1 +T*^")|r = 0), and (l+rfw = 0), are actually a superposition 
of |r = 0, Ti - 0) and \T - l,T2 -0), such that both even and odd partial waves contribute. We 
find that the former give rise to 

M^lo = -{-^/iT-fiWGs + G,)^m, (B.14) 

while the latter lead to 

M^li = -i-y/lT-fiWG,. - G,.)^m. (B.15) 

Therefore, the analog of dB.lb becomes 

X {'G.,(q' --,q--+l;-e-ko, -k) + /G,(q' - -, -q -H - - 1; -e - ^o, -k)} 

i i i 

-vl + iT] (M^v + If -Ml + ij] {M„v + kf- - Ml + irj ^ 

X {iGs{q' --,q--+l;-e-ko, -k) + /G,(q' - -, -q + - - 1; -e - fco, -k)} 

/ / / 

^ -vl + irj (M„v + If -Ml + iT] (M^v + k)^ - Ml + irj ^ 
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where again the overall factor of 2 accounts for the time-reversed diagram. Performing the P and 
k^ integrations yields 

r^wM ^ , r d\ r dH fd^qd^q' ... , 1 



: {2G.(q' - | 



k k , k^ X 
q +l;-e ,-k 



+ 2G.(q'-|-q.^-l;-.-^,-k)}T(q). (B.17) 

To generalize (IB. 5b we must replace the isoscalar two-body contribution to the jr^d scattering 
length 2T^ by the double-scattering analog 

-MT-f -f--^»i'f(p-q)->P(p), (B.18) 

J i^ny q^ 

such that the infrared-safe amplitude is given by 



^iRsafe -^^(^ ) j (2;r)3 j (2;r)3 k^P J (2;r)3 

x^W){2G.,(q'-|q-^.l;-.-^,-k) 



k^ 

2M;, 
k k , k^ . 4(27r)^^(q')*I"(q - 1) 



. 2G.(q- - ;, -q . ^ - 1; -. - ^, -k) - ^^^^J^ '\^^ " '' >F(q), (B.19) 
which finally proves ( 15.26b . 

Appendix C. Photon diagrams in heavy-pion EFT 

Appendix C.l. Electromagnetic radius of the deuteron 

As the correction due to the electromagnetic radius of the deuteron is already included in the 
experimental value for the level shift, we should subtract the corresponding NREFT contribution 
as well (in addition to ( IB. 51 )) 



.y^radius ^ 2 ri-L/2r+ ^ '- {-ie)(ie)U - -^\rl)\ 

'■* J {Inf -ko - k2/2md + ir] ka - B^t - ^-^IIM^ + irf^ '^ \^^ 6 ^ '"' 

(C.l) 
where B^i - a^poll includes the atomic binding energy and (r^) = (Sy^)"' for asymptotic wave 
functions. Performing the k^ integral yields 

J (27r)-^k- -Bat-k2/2/Uo^ 6 " ' 



y^mdius ^ ^2e'2T 



In fact, this expression is linearly divergent: the theory is only valid in the hadronic-atom regime. 
In dimensional regularization we may drop a scaleless integral, which leads to 

^mdius ^ £^^^2^2;zB yl2il^,, (C.3) 
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and subtracting this contribution from (15.34b gives 



a 



FF, radius 



= --o'(l+log4+— -)— -«\ (C.4) 



4r^d^ r e^ 



Therefore, the additional shift due to the deuteron charge radius ~ aM„/4'/^d - 0.005 is certainly 
negligible. 

Appendix C.2. Double scattering with photon exchange 
We start from (I5.27l i in the form 



......no. , 32nM^,a-f r d^k r d^l 1 I 

^d J (2;r)3J (2;r)3k2l2^" 



(rf9)+(^io) ^ -^^^' "v^p" ^ ( " »■ j " ■ i , f (k)f (k - 21) - F(2l) 

-,-\2 



ISen'-aM^i^pa-y f d^k y / „ , 4y |k 



?d 



where 

^3i j3 



r-^X'-'-ii -'"")■ <«> 



^^'"^'^ J (2;r)3(27r)3l2(p-k/2 + l)2+y2p2+y2- ^^'^^ 

Introducing Feynman parameters in the standard way, this becomes 

1 r' dx r^-' r dH ( y / y \i^ \- x . ,1-^/2 

^^C^^ r "y [^ I (C7) 

167T X vr^ Jo (j(i - 3')(i -x) + x)3/2 j (my (12 + /^(i _ h)^ + ^2)3/2 ' V ■ ^ 

with 

y(l-y)(l-x) .2 1-^ 2 .r-ox 

h - , y - y , (C.8) 

y(l-y)(l-x)+x' ^ 3;(i_3;)(i_^) + ^''' 

and hence 

/(iki)-/(0)=--^ r -^ r — — -^ — -5-^ log (1+ Ml -/z)^). (C.9) 

64;r3 J^ VT^ Jo (y(l - y)(l - x) + x)^^- \ 4y2 / 

The corresponding contribution to (IC.5l l is given by 

16ffM,(V-)2 pdz z fi dx r'dylog(l + ^(i^^i^z2^ 



r°°dz z r 

I ^- arctan - ( , , 

n^e^ Jo z^ 2 Jo Vr^Jo (y(l-y)(l-x) + x)3/2 

8QfM;,(^pa-)2 [;r2 f' dx f' dy f" dz/2 z \ 

^2^2 I Y + J„ vf3^J„ (,(i-,)(i-.) + .)3/2j„ i2i--'^t-5-ij 

xiog(i. ^^ -^i;-^\ z^) 

^g-—-- 2.916, (C.lO) 
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such that we are left with 



fiicc = TT-:, — -2.916 T 7(0) -— p<— -arctan-— 

Ti^£l \2 I £l ' ' J k4 I |k| Ay 



FF TT^^l \2 / ^2 ''"'J kM|k| 4y 



= 5-^ T" 2-916+ '' m- (C.ll) 

Now, 7(0) does not converge: the reason is that the double-scattering contribution to the ji^d scat- 
tering length itself diverges in heavy -pion EFT and we need a counterterm D to renormalize this 
diagram. Demanding that the physical double-scattering contribution in the static approximation 
^static ^g reproduced, we have 

«"^ = 4{^ - I'^^^'A^)] = 4^^ - 512.V(^P«-)^/(0)}. (C.12) 

However, the same counterterm will contribute to ((ig) and (liio) as well and needs to be added 
to ( IC.llI ). Eliminating D therein in favor of 7(0) and a^"""^, 7(0) cancels and we finally ob- 
tain ( 15.36b . 

Appendix D. Subtraction of virtual-photon effects 

Appendix D.l. nN scattering lengths 

The isospin-violating corrections to the scattering lengths calculated in 12411 involve both 
contributions due to the pion mass difference ~ e^Z, and those due to virtual photons ~ e~, where 
Z = {Ml - M^„)/2e^Fl = 0.81. (The nucleon mass difference does not enter in these channels at 



third chiral order.) Retaining only the e~ part, we obtain the following correction 124 1 

9 -) 

where g] and k'l denote the e~ piece of g] and A:J, respectively. The relation between both sets 
of LECs can be established by means of their /^-functions a-/ and ;;,. It is convenient to define 
scale-independent LECs g/ and k, by 



^ = T&(*-'-t)- «-Tife(*-'-T)' 



such that 

M^\ ^,. 77,lz=n /_ . M,r 



J-r O-i\z=0 (j , , M„\ .,. TJi\z=0( M„\ 



p. I ibn'-t^y p 

which finally leads to 






Estimating the LECs as in [24] yields the result quoted in (19.12b . 
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p \ n / p P \ " ,' P ^i / \ \' ^ \ '^ / P P ^K n / p 



n p p n p 

(a,) (a^) (a^) 



rr/ ■■, TT TT/ \Tr 



fa,; 




p n p p n p p n p 

(3e) (a,) 



Figure 7: Virtual-photon corrections to the nucleon-pole diagram mn p scattering at 0(p ) 12511 . 

Appendix D.2. nNN coupling constant 

The full set of virtual-photon corrections to the nucleon-pole diagram in n^p scattering at 



third chiral order is depicted in Fig. [7] Based on the discussion of these diagrams in 12511 . the 
shift of gc due to virtual photons can be read off from the residue of the scattering amplitude at 
s - m\. To obtain an idea how large these effects are, we consider here the diagrams (ai), (02), 
and (fls), which is motivated by the expectation that it ought to be possible to absorb this subset of 
diagrams into a simple redefinition of gc- And indeed, we find that these diagrams, together with 
the pertinent contact terms and the wave-function renormalization, represent a scale-independent 
quantity: all ultraviolet divergences cancel between loops and contact terms. In a strict chiral 
expansion, we find that (a\), (02), and (05) yield (post-renormalization) a shift in g^ of 

Ag? '^'jlj^aFl 20-, 1/ ^^U 1 1 ^' \ m^^ 



max 



The last term is present because infrared divergences only cancel at threshold. To remove these 
singularities (regulated by a finite photon mass niy in the actual calculation of the diagrams), we 



use the leading, logarithmically enhanced part of the bremsstrahlung calculated in 112511 , which ef- 
fectively eliminates niy in favor of twice the detector resolution Zsmax- Numerically, this amounts 
to 

A§? 

-P = 0.07 ± 0.03 ± 0.04, (D.6) 

An 

where we take iSmax - lOMeV and the errors are, respectively, due to the LECs and a variation 
of iimax by a factor of 2. Dropping the term due to bremsstrahlung, the shift in g^/An is reduced 
to -0.01+0.03. We stress that these estimates can by no means replace a full analysis of radiative 
corrections, but we think they can be taken as indicative of the size of virtual-photon effects. 
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